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PREFACE. 


This Treatise is in part a republication of those portions of 
my work on Mechanical Philosophy which treat of Attractions, 
Laplace’s Functions, and the Figure of the Earth. 

The first edition, issued last year, consisted of a small 
number of copies. In the present issue the last Chapter ha3 
been rearranged and in part rewritten: other improvements 
have been made. 

The disappearance of the Mechanical Philosophy has re- 
moved from the student — at any rate for the present, as no 
other work has yet appeared in English to supply the 
want — one subject of great importance and high interest, 
which that work first introduced into the University ; I 
mean Laplace’s Coefficients and Functions and the calcula- 
tion of the Figure of the Earth by means of his remarkable 
analysis. The late Professor O’Brien subsequently published 
a Tract on the same subject; but it was incomplete. A 
Fourth Edition of Mr Airy’s Tracts has been recently 
published, and in these is a treatise on the Figure of the 
Earth. But he adheres by choice (as stated in his Preface) 
to the “ geometrical and quasi-geometrieal methods.” There 
is still room, therefore, for the present Treatise; as no 
student of the Higher Branches of Physical Astronomy 
should be ignorant of Laplace’s Analysis and its results — 
“ a calculus,” to use Mr Airy’s language,. “ the most singular 
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PREFACE. 


in its nature, and the most powerful in its application that 
has ever appeared*.” 

There are problems in the Figure of the Earth which the 
geometric and quasi-geometric methods cannot touch, and of 
which the student must remain ignorant, if he is ignorant 
of the method of potentials. 

It has been my endeavour to put the well-known difficulty 
in Laplace’s analysis, arising from the use of a discontinuous 
function, in the clearest light, that the student may understand 

both what it is and how it is overcome. I have made use of 

« 

Professor Stokes’s valuable Paper in the Cambridge Philo- 
sophical Transactions of 1849 on the “Variation of Gravity 
at the surface of the Earth.” I have also introduced some 
Propositions on the Geodetic Method of determining the 
Figure of the Earth, suggested by an acquaintance with 
the circumstances of the Great Trigonometrical Survey of 
India, and by the volume of the Ordnance Survey of Great 
Britain and Ireland recently published. 


Calcutta, i86j. 


J. H. P. 


* See Article on Figure of the Earth, in the Encyclopedia Metropolitana, 
p. 19s. 
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ATTRACTIONS AND LAPLACE’S 
FUNCTIONS. 


1. The Law of Universal Gravitation teaches us, that 
every particle of matter in the universe attracts every other 
particle of matter with a force varying directly as the mass of 
the attracting particle and inversely as the square of the distance 
between the attracted and the attracting particles. Taking 
this law as our basis of calculation, we shall investigate the 
amount of attraction exerted by spherical, spheroidal, and 
irregular nearly-spherical masses upon a particle, and apply 
our results in the second part of this Treatise to discover the 
Figure of the Earth. We shall also show how the attraction 
of irregular masses lying at the surface of the Earth may be 
estimated, in order afterwards to ascertain whether the irregu- 
larities of mountain-land and the ocean can have any effect on 
the calculation of this figure. 

CHAPTER I. 

ON THE ATTRACTION OF SPHERICAL AND SPHEROIDAL 

BODIES. 

Prop. To find the resultant attraction of an assemblage of 
particles constituting a homogeneous spherical shell of very smaU 
thickness upon a particle outside the shell: the law of attraction 
of the particles being that of the inverse square. 

2. Let 0 be the centre 
of the shell, P any particle 
of it, OP= r , dr the thick- 
ness, C the attracted particle, 

/ POC = 6 ; mPMn a plane 
perpendicular to OC, <f> the 
angle which the plane POG 
makes with the plane of the 
paper, PC = y. 

p. A. * 1 



Digitized by Google 



2 


ATTRACTIONS. 


The attraction of the whole shell evidently acts in CO. 

Let OP revolve about 0 through a small angle dd in the 
plane MOP’, then rdd is the space described by P. Again, 
let OPM revolve about OC through a small angle d<f> , then 
r sin d d<f> is the space described by P. And the thickness of 
the shell is dr. Hence the volume of the elementary portion 
of the shell thus formed at P equals rdd . r sin dd<f> . dr ulti- 
mately, since its sides are ultimately at right angles to each 
other. 

Then, if the unit of attraction be so chosen, that it equals 
the attraction of the unit of mass at the unit of distance, 
the attraction of the elementary mass at P on C in the 
direction CP 


_ pr_ sin djJrddd<f> * ^ ^ density of the shell ; 

... attraction of P on 0 in CO = £ c - r co, fl 

/ y 

We shall eliminate 6 from this equation by means of 

y* = c s + r* — 2cr cos 6, 

. a dd y Q w* + c s — r s 

dry cr 2c 

.’. attraction of P on C in CO — ^1 + - dy d<f). 

To obtain the attraction of all the particles of the shell we 
integrate this with respect to <f> and y, the limits of <p being 
0 and 27 r, those of y being c — r and c + r ; 

attraction of shell on C = J J ^1 + ^~ r ^ dydcf> 


irprdr 


£C(" £ r : )*- s £p'+* 


_ iirpPdr _ mass of shell 
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SPHERICAL SHELL ON INTERNAL POINT. 


3 


This result shows that the shell attracts the particle at C 
in the same manner as if the mass of the shell were condensed 
into its centre. 

3. It follows also that a sphere, which is either homo- 
geneous or consists of concentric spherical shells of uniform 
density, will attract the particle C in the same manner as if 
the whole mass were collected at its centre. 

Prop. To find the attraction of a homogeneous spherical 
shell of small thickness on a particle situated within it. 

4. We must proceed as in the last Proposition; hut the 
limits of y are in this case r — c and r + c ; hence, 

attraction of shell = — p -f r ~ f f 1 — t dy ' 
c J T-c \ y / 

= ier*:( 2 c _2c) = 0; 

c 

therefore the particle within the shell is equally attracted in 
every direction. 

5. This result may easily he arrived at geometrically in 
the following manner. Through the attracted point suppose 
an elementary double cone to he drawn, cutting the shell in 
two places. The inclinations of the elementary portions of 
the shell, thus cut out, to the axis of the cone will he the 
same, the thickness the same, hut the other two dimensions 
of the elements will each vary as the distance from the at- 
tracted point; and therefore the masses of the two opposite 
elements of the shell will vary directly as the square of the 
distance from that point, and consequently their attractions 
will be exactly equal, and being in opposite directions will 
not affect the point. The whole shell may he thus divided 
into pairs of equal attracting elements and in opposite direc- 
tions, and therefore the whole shell has no effect in drawing 
the point in any one direction more than in another. 

6. The results of these two Propositions are so simple 
and beautiful, that it is interesting to enquire whether these 

1—2 
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properties belong exclusively or not to the law of the inverse 
square of the distance. To determine this is the object of the 
four following Propositions. 

Prop. To find the attraction of a homogeneous spherical 
shell on a particle without it; the law of attraction being repre- 
sented by <f>(y), y being the distance. 


7. The calculation is exactly analogous to that given 
above : we have only to alter the law of attraction. Then 
attraction on C in CO 




+ c* — r s ) <f> (y) dy (integrated by parts) 


= ^ IV+C-C)H> (») <&-*/&#(»)<&} ijf] 
(/ 


= 7rpr f- {(y* + d — r*) <f> t (y) - 2 i|r (y) + const.} suppose, 
c 

(c+r 1 c v 1 I 

=2t rprdr ft (c+r) - ^ ft (c+r) - — ft (c-r) + ^ ft (c-r) j- 

=2irprdr * ft (■> + 0 - + (c - r)) _ 

this latter form being introduced merely as an analytical 
artifice to simplify the expression. 

Prop. To find the attraction of the shell on an internal 
particle, with the same law. 


8. The calculation is the same as in the last Article, 
except that the limits of y are r — c and r + c\ 


.*. attraction = 2-irprdr 


7 ^<f> l (r+c)-^(r + c) 


+ ~ ft (*■ “ c) + ^ ft (r - c)j 


= ^dc\ 


d (yfr (r + c) — yfr (r — c) 


} 
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Prop. To find what laws of attraction allow us to suppose 
a spherical shell condensed into its centre when attracting an 
external point. 

9. Let <f> (r) be the law of force ; then if c be the distance 
of the centre of the shell from the attracted point and r the 
radius of the shell, and 

'lr(r) = f{rjcf>(f) dr}dr, 
then the attraction of the shell 




But if the shell be condensed into its centre, the attraction 
= 47rpr 2 dr <f> (c) ; 




(y/r (c + r) — yfr (c — r) 


} 


_ d /difrc r d 3 yjrc r 3 1 \ 

dc\dc c + dc 3 c 1.2.3 + "V 

__ , . . d (d 3 ^rc r 3 1 \ 

= 2r^(c) + 2^(^- T - 0+ ...j ; 

.% ^ d = 0, whatever r be; 

• — (\ d ^ c \ - o £ a _ o 

• • dc\c dc 3 )- ’ dc \e dd ) ~ ’ *•* 

But ^ = cf<f> (c) dc, = J(f> (c) dc + c<f) (c), 


dc 3 ~ 2< P C + C dc ’ 


therefore by the first of the above equations of condition 

2 . d<bc . „ . 

- <pc + -f- = const. = 3 A, 
c T dc 
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and. multiplying by c* and integrating 
c* <f> (c) = Ac 3 + B, 

A and B being independent of c, 

<t> (c) = Ac + ^ . 

This is the most general solution of the first of the equations 
of condition for -»|r (c), and it satisfies all the rest. Hence the 
only laws of attraction which have the property in question 
are those of the direct distance, the inverse square, and a law 
compounded of these. 

Prop. To find for what laws the shell attracts an internal 
point equally in every direction. 


10. When this is the case 


d {y]r (r + c) — yfr (r — c)| _ A 
do \ c J ~ U ’ 

+ ... = — A, 


d^rr (Pyjrr c 4 
dr dr 3 1.2.3 


whatever c is, A being a constant independent of c ; 

• QH- A n 

** dr ~ dr 3 ~ U ’ "• 


These conditions are all satisfied if the first is : this gives 
rJ<f>(r)dr = -A, <}>(r) = ~, 

and therefore the inverse square is the only law which pos- 
sesses this property. 


11. The form of the Earth and of the other bodies of the 
Solar System differing from the spherical, and more resem- 
bling the spheroidal, it is desirable to find the attraction of a 
spheroid upon an external and an internal point. 

Prop. To find the attraction of a homogeneous oblate sphe- 
roid upon a particle within its mass; the law of attraction being 
that of the inverse square of the distance. 
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12. Let a, c be the semi-axes; the minor axis 2c coin- 
ciding with the axis of z : then the equation to the spheroid 
from the centre is 


* ±y 



Let fgh be the co-ordinates to the attracted particle, which 
we shall take as the origin of polar co-ordinates, j 

r — radius vector of any particle of the attracting mass, 

0 = angle which r makes with a line parallel to z, 

(f> = angle which the plane in which 0 is measured makes 
with the plane xz ; 

x=f+r sin 0 cos(f>, y = g + r sin 0 sin <f>, z = h + r cos 0, 
and the equation to the spheroid becomes 

(f+r sin 0 cos <f>)*+ (g + r sin 0 sin <£) 8 _ (/i + r cos#) 8 t 
_ — H „ — 1, 


or r 


(^+=?-Ve 


+ 2 r ® C0S ^ +.7 sin $ sin ^ cos 
d 2 c' 




put 


1 f+lf ff. 

a* c 8 ’ 

sin 8 0 cos 8 0 „ 

— r- + = K, 

a c 


f sin 0 cos (f) + g sin 0 sin tp t k cos 0 ^ 

-- T o > 

a c 

and F*+K(l-£^-^=ff, 

then K*r t +2KFr + F* = H, 

and the values of r are 

-F+^H) - F -V(H) 

r - K and r ^ 
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The volume of the attracting element = r 3 sin 9drd0d<f> as 
in Art. 2 : let p be the density of the spheroid. Then the 
attraction of this element on the attracted particle is 

p sin 0drd9d<f> : 

and the resolved parts of this parallel to the axes of xyz are 
p sin 3 6 cos (f>drddd<f>, p sin 3 6 sin <f)dr d9d<f>, 
p sin 6 cos ddrddcfy. 

Let A, B, C be the attractions of the whole spheroid in the 
directions of the axes, estimated positive towards the centre 
of the spheroid. Then these etpial the integrals of the attrac- 
tions of the element; the limits of r being — r and r", of 6 
being 0 and v, of <f> being 0 and tt. Hence 


A = — f f f p sin 3 6 cos <f)drd0 d<f>, 

B = — f I f p sin 3 0 sin <j)drd0d<]), 

J —rJ o^O 

C — —j j f p sin 0 cos 0 dr d0d<f>, 

A = —p( f (r" + r') sin t 0 cos<f>d0d<f> 

J OJ 0 

n w P 

-^sin 3 9 cos <j>d9d<f>. 


Now it is easily seen that if It (sin a, cos 3 a) be a rational 
function of sin a and cos* a, then 

f R (sin a, cos 3 a) cos a d% = 0. 

J 0 

Therefore by substituting for F and K we have 


A= 



sin 3 9 cos 3 <f> d9d<f> 
c 3 sin 3 9 + a 3 cos 3 9 
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HOMOGENEOUS SPHEROID ON INTERNAL POINT. 9 
„ f n sin 8 Odd r s [* (1 - COB*0) sin 6d6 ^ 

= 7r > J 0 + ~ W* J 0 c 2 +(a a - c ! ) cos 2 0 

= 7r/p- j « , (* Sin >f on „' a ~ si p 4 d9 

^ {- 777 =?' tor ‘ (rr 1, ”4 +m ‘ e + ' co ” st ) 

<£—<? (cVa 2 — c 2 C J o 

„ ^ fVl-e'\ e 1-e 2 ) 

= 2 ’*{^- t “ vTr?“— } 

„ ^ 1 - e 2 } 

= 2ti/p |— J— sin e - — j . 

In the same manner we should find that 

fVTT? 1 - e 2 ) 

B=2irgp\ — j— sin e ^r-j. 


Also C — 2pJ j ~ sin 0 cos 6d6d<f> 

n i a [ n [ ,r sin 6 cos‘0d0dd> 

= n‘O + a-cos'9 

„ 7 «* [’ • d c 2 sin 0 1 m 

= 27rpA -3 s Ism 0 — a , / 8 5-5 ' 

^ a — cj 0 c 2 + (a 2 - c 2 ) cos 2 0 J 

. 7 a* (, c 4 .jVa 2 — c 2 [ 

r o 2 -c 2 l Va 2 - c 2 C ) 

, 2 fl Vl-e 2 . 1 

= 47rpA|- s ^ — sin ej. 

* If the spheroid be prolate, c is > a and the denominator of this must be 
written c 9 — (c* — a*) cos* 6, and the integral would involve logarithms instead 
of circular arcs. 
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13. We gather from these expressions, that the attraction 
is independent of the magnitude of the spheroid, and depends 
solely upon its ellipticity. Hence the attraction of the sphe- 
roid similar to the given one, and passing through the attracted 
particle, is the same as that of any other similar concentric 
spheroid comprising the attracted particle in its mass. Hence 
a spheroidal shell, the surfaces of which are similar and con- 
centric, attracts a point within it equally in all directions. 
This property can be proved geometrically exactly as in 
Art. 5. 

14. If we put the ellipticity of the spheroid = e, and sup- 
pose e so small that we may neglect its square, we have 



If we had taken an ellipsoid instead of a spheroid, the ex- 
pressions would not have been capable of integration. 

15. If we had attempted to find the attraction on an exter- 
nal particle according to the process of the last Article, we 
should have fallen upon expressions which no known methods 
have yet integrated : and therefore we are unable by any di- 
rect means to obtain the attraction of a spheroid on an external 
particle. Mr Ivory has, however, devised an indirect method 
of obtaining it, which we shall now proceed to develop. He 
has discovered a theorem by which the attraction of an ellip- 
soid upon an external particle is shown to be proportional to 
that of another ellipsoid, dependent on the first for form and 
dimensions, upon a particle internal to it, and therefore (in the 
case of a spheroid, or ellipsoid of revolution) determinable by 
the last Proposition. 

Prop. To enunciate and prove Ivory's Theorem. 

16. Let ^ + |j + - a = 1 > a nd^i + ^ + y = 1, 




SPHEHOID ON EXTERNAL POINT. IVORY’S THEOREM. 11 


be the equations to the surfaces of two ellipsoids having the 
same centre and foci : then 


a 2 -b 2 = a 2 -p ! , a 2 — c* = a 2 — 7* (1). 

Let fgh, f'g'h' be the co-ordinates to two particles so situ- 
ated on the surfaces of these ellipsoids that 


/ 

f 


a 


q _ b h 

g'-p’ r 


0 

V 


.( 2 ). 


Also since (fgh) (f'g'h') are points in the surfaces of the first 
and second ellipsoids respectively, we have 


f f h 2 , f' 2 g' 2 h' 2 , /0 , 

~ 2+ b 2+ c 2=h a 2+ W + f =1 (3) * 


Then the attraction of the first ellipsoid parallel to the axis of 
x on the particle at the point (f'g'h') on the surface of the second, 
is to the attraction of the second ellipsoid on the particle at the 
point (fgh) on the surface of the first in the same direction, as 
ah : a/9, the law of attraction being any function of the dis- 
tance: and similarly with respect to the axes of y and z. This 
is Ivory’s Theorem. 

We shall, for convenience, represent the law of attraction 
by the function r<f>(r 2 ), r being the distance. 

The attraction of the first ellipsoid on the particle (f'g'h') 
parallel to the axis of z 

= pfff(h' -z)<f> {(/' - x) 2 + (g -y) 2 + (h! - z)’} dxdydz, 


the limits of z are — c 
the limits of y are — b 





and the limits of x are — a and a 


=pim(f-xy+(g'-yy+(h'+z) 2 } 

-ir{(f'- *)* + (f - y)*+ (*' - *)*}] fady 
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12 ATTRACTIONS, 

between the specified limits : 

it must be remembered that in this expression 

-V( l -S-0- 

but we do not substitute this value merely that the function 
may be preserved under as simple a form as possible. Now 
put x=ar, y = bs, z = ct, then the attraction 

= P<dfS {(/' - ar)' + {g - bs)' + (Ji - c<)*} 

— ^ — ar Y + iff' ~ hs)* + [h' + c<) a }] drds, 

the limits of s being — V( 1 — r') and V(1 — r*), and those of r 
being — 1 and 1 : also t = *J(\ — r* — s *). 

Now (/' - ary + {g - bs)' + {K ± ct)' 

=f" + g" + h! 2 — 2 ( far + g'bs + h'ct) + oV + b's' + cV, 

substituting for h' 2 by (3) and for f, 

=/'* ( 1 - ^1) + /* ( ! + •/-» (f'ar+g'bs ± h'ct) 

+ (a'-c')r'+(b'-c r )s , + d, 
eliminating f'g'K by (2) and making use of (1), 

= ^ (a* - c*) + 1 5 {b'-c')+<?~ 2 (/ar +y/9s i A7O 

+ (a 2 — 7 s ) »•* + (^* — 7 s ) s a + 7* 

=/* + / + h* - 2 (/ar + gj3s + hyt) + aV + + 7V, by (3) , 

= if- ar )' +(g- fis)' +{h± 7 1). 

Hence the attraction of the First Ellipsoid on ( fgh ') 
parallel to z, 

= P ah // ( if- ar Y + iff- @ S Y +(h + 70’} 

- ir [if~ °- r Y + (9 ~ P S Y + (* - 70 s }] drds 
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= —x attraction of Second Ellipsoid on (fgh) in the same 
direction. 

The same may be proved for the attractions parallel to the 
other axes: and consequently the Theorem, as enunciated, 
is true. 


We may observe that one of these ellipsoids must neces- 
sarily be wholly within the other. For if not, the points in 
which they cut each other lie in the line of which the equa- 
tions are 


-, + ^+3 = 1 and^ + ^+i = 

a b c a p 7 


1. 


Suppose a less than a; the points of intersection must 
satisfy the equation 


x 




and this by (1) becomes 


( 0 + (£)-* 


an equation which can be satisfied only by x = 0, y — 0, 
z — 0. But these do not satisfy the equations above; and 
therefore the surfaces do not intersect in any point. 

To find the attraction of any ellipsoid of which the semi- 
axes are a, b, c upon an external point ( fgh ') by the help of 
this Theorem, we must first calculate the attraction of an ellip- 
soid of which the semi-axes are a/3y, determined by equations 
(1) and the second of (3), on an internal point ( fgh), f, g and 
h being given by equations (2). And then tlie attractions 
required will be those multiplied by 


be clc cib , • i 

respectoely - 
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LAPLACE’S COEFFICIENTS AND FUNCTIONS. 


17. In the present Chapter we shall develop the properties 
of those remarkable quantities which have received the name 
of their great discoverer, under the designation of Laplace’s 
Coefficients and Functions. To do this it will be neces- 
sary to anticipate the subject of the following Chapter, and 
to bring in here a Proposition which should properly stand at 
the head of that division of this treatise. 


Prop. To obtain formulae for the calculation of the attraction 
of a heterogeneous mass upon any particle. 

18. Let p be the density of the body at the point (xyz) ; 
fgh the co-ordinates of the attracted particle ; ana, as before, 
suppose that A, B, C are the attractions parallel to the 
axes x, y, z. Then 

p ( f — x) dx dy dz 


" ///{(/- xy + {g- yY + (h - zf p ’ 
iff p (<7 — y) dxdudz 

"iJi 

P- f f f P (* ~ z ) dxdydz 

0 "iii{(/-^+(y-y) i +(A- z rr 


B 


the limits being determined by the equation to the surface of 
the body. 


Let 



pdxdydz 

{(/-*)•+(*->+ (*■ 



c =- 


dV 
dh • 


19. It follows, then, that the calculation of the attractions 
A, B, C depends upon that of V. This function cannot be 
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calculated except when expanded into a series. It is a function 
of great importance in Physics : and, for the sake of a name, 
has been denominated the Potential of the attracting mass, as 
upon its value the amount of the attractive force of the body 
depends. 


20. As the axes and origin of co-ordinates in the previous 
Article are altogether arbitrary, it follows that if r be the 
distance of the attracted point from any fixed point in the 
attracting body, then the attraction in the line of r, towards 

. . , dV 

the origin of r, = — - 3 - . 


d*V ddV 

Prop. To prove that -jj? + - 


d'V 
' did = 


= 0, or — i7rp, ac- 


cording as the attracted particle is not or is part of the mass 
itself ; p being the density of the attracted particle in the 
latter case. 


21. By differentiating V, we have 
dV iff —p (f—x) dxdydz 

k f-xr+ig-yr+ih-zy]^ 

<PV _ fff p {2 (/- x)' - (g —yf— {h-zf) dxdydz 
df~)J) {(f-xy+(g-y)'+{h-zY}l 

In the same manner we shall have 

d' i V_ fffp ( 2 (. 9 - y? - if - x )* - ( 7i - 2 )'l dxdydz 

df-JJI {{f-xy+(g-yy + (h-zf}i 

d‘V_ ff[ p{2(h- »y - {f- xy - {g - y) 8 } dxdydz . 
dh'-JJJ { (f-xy+(g-yy+(h-zy}i 

d'V d'V d'V fff Ox dxdydz 

dr + W + ~did - J 1 J {(/- xy + (g - y y + (h - «)*]! * 

When the attracted particle is not a portion of the attract- 
ing mass itself, then xyz will never equal fgh respectively, 
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and consequently the expression under the signs of integration 
vanishes for every particle of the mass : 


d'V fV d'V_ 
*** df + dg* + dli* ~ 


This equation was first given by Laplace: and Poisson was 
the first who showed that it is not true when the attracted 
particle is part of the attracting mass. In that case the deno- 
minator of the fraction under the signs of integration vanishes, 


and the fraction becomes - , when x =f y = g, z = h. 


To determine the value 


, d'V t d'V a d'V. . 
of - 7 s + -77 + -jti in that case, 
dj dg dh 


suppose a sphere described in the body, so that it shall include 
the attracted particle ; and let V — U+ U ' , U referring to the 
sphere, and U' to the excess of the body over the sphere. 
Then, by what is already proved, 


d'U' , d?U' , d*U' ^ 
df + dg 1 + dh* 5 


<PV d 2 V d*V L d'U d'U d'U 
*’• df + dg* + (Th ” df* + dg* + dh* ' 

The centre of the sphere may be chosen as near the 
attracted particle as we please ; and therefore the radius of 
the sphere may be taken so small that its density may be 
considered uniform and equal to that at the point ( fgh ), which* 
we shall call p. 

Let fgh' be the co-ordinates to the centre of the sphere ; 
then the attractions of the sphere on the attracted point 
parallel to the axes are, by Art. 3, 


47 rp 

~ 3 ~ 


(/-/'), 


-3- (£-£'). 


A-TTp 


(h-h') f 


or 


dU 

df' 


clU 
dg * 


dU 

dh' 


by Art. 20. 
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d'U , d*U d'U . , 
*’* df a + dg t + dh a ~ i7rp ’ 

d*V d*V cPV A , 
" df + dff + dh a ~ i7rp ’ 


when the attracted particle is within the attracting mass. 


22. It may be shown by precisely the same process as in 
the previous Article, that 

d a R (PR d‘R_ 

df a+ dg a + dh a ~ 

where R = {{f-x) a + (g-y) a + (h- z) a p, 

the reciprocal of the distance of any point of the body from 
the attracted particle. 

Prop. To transform the partial differential equation in R 
into polar co-ordinates. 


23. Let rOco be the co-ordinates of {fgh), and r'0'co of 
(xyz), the angles 0 and 0' being measured from the axis of z ; 
co and co' being the angles which the planes on which 0 and 
0' are measured make with the plane zx. Then 

f=r sin# cos G), g — r sin 0 sin co, h = r cos 0, 

x = r' Bin 0“ cos cd, g = r sin 0' sin co', h' = r cos 0'. 

These are the same as 

r'=f*+f+h*, cos 0 = ^=^=, tana, (1) ; 

dR _ dR dr dR d0 dR dm 
*’• ~df~df df + d0 cf + fo ~df' 

d*R _ d dR dr d dR d0 d dR dco 
dff df dr df+dfdOdf+dfdoodf 

dR d'r dR d*0 dR cPco 
+ dr df t+ d0 df + dio~df 
P. A. 2 

ST 
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_ d ' R dr l d * E d0> dw ' 

~ dr* df + d0* df + da? df 

d*R dr dd d*R dr da> d*R dd da> 
^ drdd df df+ dr day df ~df+ ddda> df df 


dR <Pr d_Rfd dR d?a> 

+ dr df 1 + dd df + da> df ’ 

(PR PR 

The expressions for and are of the same form. 

These three must be added together and equated to zero. 
When this is effected the formula (1) make 


the coefficient of = 


dr * 

d f 


dr * dr* _ 

+ df + W~ 1 ' 


, _ . e d*R dd 1 dd 2 dd 1 

the coefficient ot -rz= = -rs + t* + jr, = 

• dd df dg dh 

, ~ . c d*R da? da? da? 

the coefficient ot —r—. = -77* + -5-* + -yy, '■ 
da? df dg dh 


I 


Vsin’0’ 


, . . d*R dr dd dr dd dr dd 

the coefficient of ^ - 2 J^+ 2 dgf} + 2 dhTh~°' 

the coefficient of = 2 ^^ + 2 ^-^ + 2^^=0, 
drda) df df dg dg dh dh 

the coefficient of — -2 — — + 2 — d< ° + 2— — -0 
the coefficient ot ~ 2 df df +2 dg dg + 2 dh dh ~ °* 

,, ~ . , c dR d*r d?r d?r 2 

the coefficient ot r = + -n + -m = - » 

dr df dg dh r 


the coefficient of + ZJT* = 


d % d 

df 


d?d 

dg* 


d*d 

dh* 


cos d 
Tsm 0 ’ 


i «... ?dR da) d to (Ta> 

the coefficient ot -j- = ~rri + -y-y + = 0. 

da) df dg dh 


d*r 
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Hence the equation in R becomes 

(PR 2 dR 1 d*R cos 0 dR 1 <PR _ 
dr* + r dr + P d0* r s sin 0 dO r* sin*# da? ~ ’ 


d*.rR d'R cos 6 dR 1 d*R 
dr* dd* sin 6 d0 + sin* 6 da? 


= 0 . 


Put cos 6 = n, then 

d*.rR d {. dR\ 1 cPR 

r sx + i' 1 "■ * > M + a? - »• 

Prop. To explain the method of expanding R in a series. 

24. The expression for R becomes, when the polar co-ordi- 
nates are substituted, 

[r* + r* — 2rr' [pp + V 1 — g? Vl — p '* cos (to — w')}]"*, 
and this may be expanded into either of the series 


P - +P — + 

o r 1 r 

n 1 r ' 

01 F oZ + F r 3 +' 


■ "b Pf itn H* • • • 


• "b Pi "b ••• 


•( 1 ), 


where P 0 , P lt ...P ( ... are all determinate rational and entire 
functions of p, 

vr: p cos (o, and Vl — /V sin w ; 
and the same functions of p, 

Vl —gP cos a>, and Vl — p* sin (o'. 

The general coefficient P t is of i dimensions in p, 

V 1 — p* cos a>, and Vl — p* sin co. 

The greatest value of P i (disregarding its sign) is unity. For 
if we put 

pp' + Vl — p 1 Vl — p* cos (to — a/) = cos <f> = ^ ^ , 

2-2 
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then P, = coefficient of c‘ in 

(1 + c* — 2c cos <£)“*, or (1 — cz)~‘ ^1 — 
= coefficient of c‘ in 


, 1 1 . 3 , 2 . W1 , 1« , 1.8cP, 

(1 + 5 CZ + 2 T 4 cz + ••') (1 + 9 « + 


2 z 2.4r 




= 2A cos i<f> + 2B cos (i — 2) £ + ... 


u4, B... being all positive and finite. The greatest value of 
this is, when <f> — 0. Hence P, is greatest when <f> = 0. 

But then P 4 = coefficient of c‘ in (1 + c s — 2c) "* or (1 — c) 1 
= coefficient of c* in 1 + c + c* + ... + c* + ... 


= 1 . 


Hence 1 is the greatest value of P ( . It follows that the 
first or second of series (1) will be convergent according as r 
is less than or greater than r. 

To obtain equations for calculating the coefficients P 0 , P, , 
...P ( ... substitute either of the series (1) in the differential 
equation in B in the last article, and equate the coefficients of 
the several powers of r to zero. The general term gives the 
following equation : 



T^7W + ’(’ + 1 ) i>,=0, 


by integrating which P 4 should be determined*. The series 
for R would then be known. 


25. The functions P 0 , P t ...P< ...possess some remarkable 
properties which were discovered by Laplace. They are there- 
fore called, after him, Laplaces Coefficients , of the orders 0, 
1 ... % ... It will be observed that these quantities are definite 

* For the direct integration of this equation, see two Papers in the Philo- 
tophical Transactions for 1841 and 1857, by MrHargreave and Professor Donkin 
respectively. 
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and have no arbitrary constants in them. Laplace’s Co- 
efficients are therefore certain definite expressions involving 
only numerical quantities with p, and a>, p and o>'. Any other 
expressions which may satisfy the partial differential equation 
in P„ which is called Laplace’s Equation, .may be designated 
Laplace's Functions to distinguish them from the “ Coeffi- 
cients.” The fundamental properties of these Coefficients and 
Functions we shall now proceed to demonstrate. 

Prop. To prove that if Q t and F t be two Laplace! s Co- 
efficients or Functions , then I I Q t R t -dpdco = 0, when i and 

J 0 

% are different integers. 


26. By Laplace’s Equation in the last Article 

n 2lr 

QJdfdpdw 

= - [| {(> - '*■> f } + R ‘ d ^- 

By a double integration by parts 

••• Li {<> - « f j *+-Li {<• f 1 «*■ 


Digitized by Google 



22 


ATTRACTIONS. 


since when as = 0 and 27 r, each of the functions Q t , R,-, 
dR t 


dQ t 

das ’ 


has the same values, because they are functions of ft, 
vT — fi* cos as and Vl — /P sin os. 

Hence, J J Q t R t - dp das 


*' (t'+l) 

tjt + l) 


ff 

J -l J o 


Q t R{ dfi das, 


by Laplace’s Equation. 

Hence, I QiR^dfxdos = 0, when i and t" are unequal. 
J —l*' 0 

When they are the same the equation becomes an identical 
one, and therefore gives no result. 

This property is true also when t=0, as may easily be 
shown by going through the process of the last Proposition, 
Q t being Q 0 or a constant. 


Prop. To prove that a function of fi, V 1 — cos as, and 
Vl — fi* sinus, as F{ji, as), can be expanded in a series of 
Laplace's Functions; provided that F\p, as) do not become in- 
finite between the limits — 1 and 1 of fi, and 0 and 27 r of as. 

27. This very important Proposition will occupy the 
present and five following Articles. 

Let p(i + Vl — /** Vl — fi 1 cos (o) — as) =p ; then by Art. 24, 

(1 + c* —2 cp) _l = 1 + P,c + P s c* + + Pfi* + 

c being any quantity not greater than unity. 

Differentiate w'ith respect to c, 


p—c 

(1 + c* — 2 cp) 1 


-P 1 + 2P e c + 


+ iP i c , - 1 + 
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Multiply this by 2c and add it to the former equation. 

•’* (i +c 3 - 2cp)* = 1 + 3-P i c + 5 -^ > c * + + ( 2 * + J ) P&i + — 

Now c being quite arbitrary we may put it = l. Then the 
fraction on the left-hand side of this equation vanishes, except 
when p = 1 ; in which case the fraction on the left hand be- 
comes apparently indeterminate : but it is in reality infinite. 

For whenj-1. ^ a - c ^ )| = infinity, when c = l. 

When p = 1 , then 


cos («'- ®) = = . / l-2pp+p a p* 

V(i - aO (1-/0 • 

and that this may not be greater than unity we must take 
fj? + p' 2 not greater than 2pp, or (/a ~ /t') s not greater than 
zero. Hence p' = p, and therefore cos (of — o>) = 1, and <a'= o>. 


These, then, are the values of p and a> which make p = 1 . 

Hence, the series 1 + 3 P t + 5P S + + {2i + 1 ) P< + 

vanishes for all values of p and to, p and ©', except when 
p — p and a = w', in which case the sum of its terms suddenly 
changes from zero to infinity. 


28. Upon this series depends the important property of 
Laplace’s F unctions which we are now demonstrating, and 
which gives them so great a value in the higher branches of 
analysis. In consequence of the discontinuity above pointed 
out, and also because the series becomes infinite in one stage 
of the variations of its variables, it has been considered by 
some to be unsatisfactory to deduce any properties from it. 
But the latter objection is entirely removed by the fact, that 
we do not use the series in its present form, but after being 
multiplied by small infinitesimal quantities which render the 
aggregate of its terms finite, preventing their accumulating 
to an infinite amount. With regard to the objection of dis- 
continuity, there appears to be no sufficient ground for it. 
There is no question, that the property deduced (as enunciated 
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in our Proposition) is true, at any rate for rational functions 
of \t, Vl — /t* cos <», and V 1 — /x a sin o>, and is also most impor- 
tant. This objection, however, deserves to be examined with 
care, which we now propose to do in the course of our de- 
monstration. , . 


29. Multiply both sides of the last equation by the double 
element dfi'dco', and integrate, 

rllt (1— c')dfi'da>' 


ff 

J -iV 0 


(1+c - ' 1 — 2 cp) 1 


— f j {+ 3cP l +...-(-(2f+l)c i Pj+...}<f/A'dJa). 
J _x J o 


The property of Laplace’s Functions proved in Art. 26, 
shows that every term of the series on the right, except the 
first, vanishes of itself, independently of the other terms; 
and therefore (as was before intimated) the terms cannot ac- 
cumulate. The first term is 4v : and therefore the integral 
of the fraction on the left, that is 

•*'(1 -c^dfi'da,' 

o (i +c *_ 2ci ,).>- 47r - 

It is remarkable that this result is altogether independent 
of c. 



30. The truth of this may be shown also by integrating 
the fraction on the left. This cannot readily be done with 
the co-ordinates as at present chosen. But it may be done by 
a simple transformation, and a change in the way of taking 
the elements. 

Suppose a sphere of radius unity described about C the 
origin of co-ordinates. Let 6' and 
<»' he the angular co-ordinates to 
a point P, & (or cos -1 ft) measured 
from a fixed point A along a great 
circle of the sphere, and o>' the 
angle which this great circle makes 
with another and fixed great circle 
through A. Then dff . da> sin ff, 
or — du'dw, is an infinitesimal 
element of the surface of the sphere 
at P. Take D a point within the 
sphere, and let VD = c, and sup- 
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pose CD meets the sphere in Q when produced forwards, and 
in q when produced backwards. Let p and to be the co-or- 
dinates of Q. Then p (see its value, Art. 27) is the cosine of 
the angle which CP and CQ ma ke with each other : and the 
distance of P from D = Vl + c* — 2cp. Let be the angle 
which the plane CPQ makes with CAQ, that is, the angle 
AQP. By changing the origin of the angles from A to Q, 
and dividing the surface of the sphere into new elements, 
beginning from Q as the origin, the element at P, with these 
new co-ordinates cos -1 jp and will be — djpdbSf. 

By reverting to the meaning of integration we see that 
the integral under consideration = (1 — c 2 ) x limit of sum of 
all the elements of the surface of the sphere divided respec- 
tively by the cubes of their distances from D. 


But this, by the change of co-ordinates, also 


f 1 f Sir (l- c*) dpd-fr 
j_J 0 (1+c*— 2cp)*’ 


which = 2w 



d p 

(1 + C 51 — 2 cp) 1 


= 27 r 


1-c* 


1 + const.) = 2tt - — — (~ =-J— ) 

c Wl+c* — 2 cp / c \l-c 1+c/ 


= 47r, whatever value c has. This coincides with the former 
result. '' v 


31. This integration helps us to see by what process c 
disappears from the result ; and it will assist us in the latter 
part of the present demonstration. 

The quantity 1 —p is the versed-sine of the arc QP, and is 
measured along the line QCq. Let this line be divided into 
n parts each equal to h, so that nJi — the diameter =2, n being 
very large and h very small. Draw perpendiculars to the 
diameter through these divisions cutting the circle QPq in a 
series of points ; and call the distances of these points from D, 
beginning from Q, 

1 — c, s’, s", s'" fi'”- 1 ’, 1 + c. 

Suppose P is at the a:" 1 division ; then 
1 —p = x.h, 

and d (1 —p) or — dp — (x+l) h — xh = h. 
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Then by mere expansion, omitting the squares and higher 
powers of h as they vanish in the limit with reference to the 
first power, we see the truth of the following ; 


'-dp 


d(l-p) 


(1 + c l — 2cp)* {(1 — c) 4 + 2c (1 — p)) 1 

-If -- 1 1 1 

c l V(1 — c)* + 2 cxh V(1 — c)* + 2c {x + 1) h) 
c U (x) s '* 41 ) • 


By giving x its successive values from 0 to n — 1, and adding 
together all the resulting values of this expression and taking 
the limit, we have the integral with respect to p. It matters 
not in which order we effect the integration. Hence the whole 
integral 


[**[' (l-c^d^dp 

L J-A 1 + c* — 2cjp)‘ 



n being made infinitely great, 



Here it will be seen that the terms within the last brackets 
mutually destroy each other whatever be the value of c. It 
may also be observed that were this not the case, that whole 
part of the expression would vanish for the particular value 
c = 1 (which is the only case we shall have to use), wlfatever 
the value of the sum of the terms following the multiplier 
1 — c*, so long as that sum is not infinite. 
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32. Suppose now F (ji, &>’) is any function of p! and ay', 
then 


ff 


_ „j\ <*') dp! day 

1 > (1 + c* — 2cp)* 


— [ j {1 + 3-iPj + ... + (2t + 1) + •••} F (ji , ay) dfi day . 

J -1*' 0 

The reasoning above enables us now to prove that the in- 
tegral on the left-hand side = 4irF(fi, ay), which directly leads 
to the theorem we are wishing to demonstrate. 

The function F{p!, ay) at the point Q is F{p,, ay), call it F: 
let F', F" ... F w be its values at the points of junction of the 
successive elements along the great circle QPq. Then by 

multiplying the successive values of - by F, F', 

0 \8 S J 

F" ... and adding them together, we have 

F (fj. 1 , ay') dp! day 
(1 + c* — Sep) 1 ’ 

F (p!, ay') dyfrdp 
(1 + c* — 2 cp)* ’ 


fjV* 

J o 

rf*-« 

J a J -l 


l+c/ j 




n being made infinitely great. 

1 — c 1 — c 


The fractions 


— , .... diminish successively in 


s s 

value, being the ratios of QD to the successive values of DP. 
When c = 1 each of them vanishes ; and in the limit none of 
the factors F' — F, F" — F‘, ... become infinite. Hence the 

integral = I dyfr . 2 F, when c = 1 , = iirF (ji, ay) because F(p, ay) 
J 0 
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is a function of p and eo only, and is altogether independent 
of yjr. 


.*. itrF(p, a>)= f [ {l+ZP l +...^i+l)P i +...)F{p!,a)dpdcd', 
J - 1 J 0 




f *»') , 3 F(p',a>') 


47T 


p t +. 


+ 


(2i+\)F(p',co') 

47T 


p+ 



The general term of this, viz. 

zx^, *o wdo>', 

which we will call F„ is a function of p and a ; and evidently 
satisfies Laplace’s Equation in p and w, because P t does so. 
Hence, this is a Laplace’s Function, of the r th order: and the 
result is, what we were to demonstrate, that any function of 
p and co may be expanded in a series of Laplace’s Functions; 
or, 

F(ji, w)=F 0 + F l + F 2 + + F t + 


33. Those who are at all acquainted with the controversy 
which followed the first discovery of these remarkable functions 
by Laplace, will understand why we have entered so fully 
upon the subject. Laplace’s demonstration in the Mecanique 
Cileste was by no means conclusive. This Mr Ivory pointed 
out in the Philosophical Transactions for 1812 ; and in the 
Volume for 1822 he threw considerable doubt upon the 
applicability of the theorem to functions that are not rational 
and entire functions of p, Vl— p * cos <w, vnr p* sin w. Poisson 
has written much upon the subject. In the first edition of the 
author’s Mechanical Philosophy the last method of Poisson 
was followed, as given in his ThSorie MatMmatique de la 
Ghaleur; in which he effects the integration of the fraction 
on the left-hand side by the artifice of substituting for it an 
integrable, but entirely different fraction in its general form, 
but which coincides with it in the particular case for which he 
requires it in the result, viz. when c=l. In the Second 
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Edition of the Mechanical Philosophy we gave a much shorter 
proof, based upon an idea taken from Professor O’Brien’s 
Mathematical Tracts. But this also rather concealed the real 
difficulty of the case, and passed it over by an artifice. In 
the demonstration now given, we have gone to the foundation 
of the calculus, the doctrine of limits, and attempted to clear 
up all difficulty and ambiguity in the matter. 

With regard to the doubt thrown out by Ivory, alluded 
to above, it seems to be clear that theoretically every function 
can be expanded in a series of Laplace’s Functions : but if it 
be not a rational function of the co-ordinates, the number of 
terms in the series will be infinite, and if the terms be not con- 
vergent, the expansion, or rather arrangement, will be use- 
less. But this must be determined in each case. A similar 
uncertainty, requiring examination, always attends the use of 
infinite Series. 

Prop. To prove that a function of p and a> can he arranged 
in only one series of Laplace's Functions. 

34. For if possible let both these be true, 

F{p, ©) = F 0 + F t + F t + +F t + 

F(p,ay) = G 0 +G 1 +G t + + G + 

••• 0 = (F 0 -G 0 ) + ft- G x ) + + (F- G ( ) + 

and if these letters be accented when p and to' are the variables 
instead of p and &>, then 

o = (f; - g;) + (f;-g;) + + (f; - g;) + 

0= f 1 r P t (F; - G;) dp' day’, by Art. 26. 

But the principle demonstrated in the last Proposition 
shows that 

F *- a <=h\\ C (l + 3P « + - ) ( F < ~ ^') 

= f I* Pi (F; - Gi) dp' da,', by Art. 26, 
= 0, by the condition deduced above ; 
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therefore F t - O,, and the two series are terra by term iden- 
tical, and the Proposition is true. 

35. It follows from this, that if by any process we can 
expand a function in a series of quantities which satisfy 
Laplace’s Equation, that is the only series of the kind into 
which it can be expanded : and if by any other process we 
obtain what is apparently another, the terms of the two series 
must be the same, term by term, and we may put them equal 
to each other. 

36. Before concluding this Chapter, we shall explain how 
the numerical coefficients in P 0 P ... P, ... are found : and shall 
give a few examples of the truth of the last Proposition but 
one (that in Art. 32) by actual integration. 

Prop. To explain how to expand P t . 

37. By Art. 24 P t is the coefficient of c 1 in the expansion 
of the function 

[1 + c* - 2c [ftp + s/l-ft 2 Vl - ft* cos (ct> - »')}]"*, 
and is therefore a rational and entire function of ft, 

Vl — ft* cos o), and Vl — ft 2 sin co ; 
and is precisely the same function of ft, 

Vl — ft 2 cos a>', and Vl — ft ' 2 sin &>'. 

The general term of P t , viz. that involving cos n (« — «'), can 
arise solely from the powers n, n + 2, n + 4, ... of cos (a — «'). 

n 

Hence (1— ft 2 ) 2 will occur as a factor of that term: and the 
other part of its coefficient will be a factor of the form 

+ Api~* "* 4- ... + Ajt i ~ n ~ 2 ‘ +... = H n suppose. 

Hence 

« 

P t =H 0 + (1 - ft 2 ) l E j cos (oi - <»') +...+(l-/x 2 ) 2 i7 n cosn(&)-o)')-f-. . . 

If this be substituted for P t in Laplace’s Equation and the 
coefficient of cos n (a — a') be equated to zero, we obtain a eon- 
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dition from which to calculate the arbitrary constants we have 
introduced. This condition, after reduction and arrangement, 
is as follows : 

0= (<•-») (<■+ »+ 1) H. (1 — #»*)* + j(l -y.y” d M . 

Substituting in this the series which //„ represents, and 
equating the coefficient of the general term (1 — to 

zero, and reducing, we arrive at the formula 

. _ (i — n — 28 + 2) (» — n — 2s + 1) 

* “ 27(2i^s+T) A, ~ l ' 

By making s successively equal 1, 2, 3 ...we have A 1 A 1 ... 
in terms of A 0 . Let these be substituted, and we have the 
coefficient of cos n ( &> — o> = 

Ml- A f /*■-”+ } , 

call this Agfiji). The coefficient A 0 is a function of /i', but is 
independent of n : and because P ( is the same function of /a' 
that it is of /*, it follows that A 0 = a n f{ji), where a, is a nume- 
rical quantity : and the coefficient of 

cos n (to - to') = a*/(/0/ (/*)• 

To find a n we must compare the first term of the ascending 
expansion of a n f (/*')/ (p) in powers of n with the correspond- 
ing term in the coefficient of c‘ in the actual expansion of 

[1 + c* — 2c [fiii 4- Vl — fA Vl — fi* cos (to — o>')}]-*. 

This leads to the following result : 

„ (1.3.5... (2i— 1)| s i (i — 1) + 1) 

“*“1 1 . 2 . 3 . . . t j (» + l)(* + 2 ) ...(* + «) J 

this applies when n — 1, 2, 3..., but evidently not when 
n — 0 : a 0 is found by equating coefficients to be 


[ 1 . 2 . 3 ... ( 2 »- 1 ) 1 * 
| 1 . 2 . 3...t J* 
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We have now the complete value of P i in a series ; it is as 
follows : 


= ( 1.3.5... (2t- 1) ] 1 
* | 1.2.3...* ) ■ 

2(2*-l) M + ' 


2 (2* — 1) 


'h «(»-!) >i-» , « (» - 1) 

2 (2* — 1) “ 2 (2t — 1) 


(«-8) («-3) 

4 (2* -3) 

( t -_2) (t — 3) 
4 (2* — 3) 


/* <-4 — &c...j 

/P--&C...J 


+ 2 cos (© — ©') 


*' 

* + 1 


X(l-/0‘ {/.-*; 

+2 cos 2 (©—ft)') 
x(l-/*’) 

xM{r- 

+ &C....1. 


(t-l)(f-2) (»-l)(t-2) (*-3)(*-4) 

2(2*— 1) ^ 2(2t— 1) 4(2*— 3) 

(*— 1) (*~2) (t— 1 ) (*— 2) (*— 3) (*—4) 

2(2*— 1) ^ 2(2*— 1) 4(2t— 3) 

f(f-l) 

(t+l)(t+2) 

(«-2)(«~3) ~ , (<— 2) (t— 3) (*— 4) (*— 5) 

2(2t— 1) ^ 2(2t— 1) 4(2*— 3) 

(»-2)(»-3) , (t— 2) (*— 3) (*-4)(*-3) 

2(2*— 1) ^ 2(2*— 1) 4(2*— 3) 



38. The following numerical examples are written down 
for convenience of reference : 


(1) Pj = fifl + Vl — fl* Vl — /*'* cos (© — ©'). 

(2) P s =j (p’*- - (l-/**)V(l-/’)Vco3 (©-©') 

+ £ ( 1 -A**)( 1 -/*'*) C °82 (©-©')} . 
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( 3 ) = 

+ 1 (i- -rWif- i) U V ) 1 (/*’- 5 ) cos («-«■) 

+ ^(1— /.'*)^t'cos2(oj— a»')+-~ /i' 3 ) , cos3(ffl— o)') • 

&c. = &c. 

39. The following are some examples of expanding a 
function in a series of Laplace’s Functions, by an application 
of the formula 

F t . = ^ f f W F(jj,', co') Pidfi'do)', 

proved in Art. 32. 

Ex. 1. Arrange a + bfj, 2 in terms of Laplace’s Functions. 
Here F(ji, ©') = « + bfi" 2 . First put i— 0, P a = 1 ; 

F 0 = ~J J (a + bfi' 2 ) d/M da)' — (a + b/i' a ) dfi 

= ^ (a/ + ^ V* + const.) =0 + ^6. 

Again, put t = l, P x is found in the last Article. 

.*. P = ~f f (a+bfjb' 2 ) [fx.fi'+'d l—fjd'd 1— ft 2 cos ((o—o)')}dfi'd(o' 

_jJ q 

— j (a+b/j,' 2 ) {/z/z. 00'— V 1— /j? V 1 —/a 2 sin (o>— &>')} dfi , 
between the proper limits, co' = 0 and a = 2i r, 

= | /Wl f A* (f V + J v) . 

between the limits /z' — — 1 and fi = 1, =0. 

P. A. 3 
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Next, put i = 2, and substitute for P s from the last Article. 

F ' “ h 0> +M {? ('‘'"I) s) + 

+ AT cos 2 (&)— cd')| dfi’day 1 

= f (^ - 1 ) /i{- 1 “ + ( “ - 1 

4V4)B a+ ! ( “-5 J > + H =i ('‘ , -f)- 

Hence the function a+bp? stands as follows, when arranged 
in terms of Laplace’s Functions, 

and consists of two Functions, of the order 0 and 2 respectively. 
The above is a long process to arrive at this result. It might 
have been so arranged at a glance. But the calculation has 
been given as an example of the use of the formula, which in 
most cases is the only means of obtaining the desired result. 

Ex. 2. Arrange 49+30/4+3//.*+ (40+72/*) cos (a-a) 

+ 24 (1 - /**) cos 2 (m - a) in terms of Laplace’s Functions. 

The result is 50 + {30/4 + 40 Vl - cos (to - a)} 

+ {3/4* -1 + 72/4 Videos (a) - a) + 24 (1 - pF) cos 2 (a -a.)}, 
consisting of three functions of the orders 0, 1, 2. 

Ex. 3. Let the function be 

1 + \/2 — 2/4* cos (a) + 0 ) + ^ (1 — /4 s ) cos 2 (&) + a). 


Digitized by Google 



EXAMPLES. 


35 


The first term is a Laplace’s Function of the order 0, and the 
second and third terms taken together are one of the second 
order. 

Ex. 4. Let 1 — (1 — i?) cos 2 u> be the function. The arrange- 
ment is 

i + 1 ~ I) - (1 ~ ^ coa 2w }’ 

or, which is the same, 

| + ||- cos 2 o>| . 


3—2 
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ATTRACTION OP BODIES NEARLY SPHERICAL. 


40. As tlie Earth and other bodies of the Solar System 
are nearly spherical, and yet may not be precisely of the 
spheroidal form, it is found necessary in questions of Physical 
Astronomy to calculate the attraction of bodies nearly sphe- 
rical. It is in these calculations that the value of the Functions 
we have been considering in the last Chapter is seen. 


If r'6'co' be the co-ordinates to any element of the attracting 
mass, p be its density, and cos & = p!, then the mass of this 
element 


= p’dr'r'dffr sin 6' dm' = — p'r" 2 dr'dp.'dco', 


and the reciprocal of the distance being R, by Art. 18 and 24, 
the potential V 


+ ••• +p ‘£ + ■■■)*#*'■• 
r J I J p\ P * r ' + P S + P * 7 + --- + Pi ~ji=i + -••')dr'dp,'d( 0 ', 


according as r, the distance of the attracted point from the 
origin, is greater or less than r . We shall proceed soon to 
use these formulae ; but we must first find the value of V for 
a perfect sphere. 

Prop. To calculate the value of V for a homogeneous sphere. 


41. Let the centre of the sphere be the origin of the polar 
co-ordinates (r'p.'ed) to any element of its mass, and the line 
through the attracted point be that from which the angles are 


Digitized by Google 


BODIES NEARLY SPHERICAL. 


37 


measured, and p the density. Then — pr'dr'dp da>' is the 
mass of the element : its distance from the attracted point 

= Vr 2 + r' a — 2 rr cos to. 

Hence, a being the radius of the sphere, 

V 


pr' 2 dr dp! da' 

o Vr a + r' a — 2rr'p 


dr', 


= 2rp FJ[^^'~^ P !yr{^ +r '’- 2 " y 

from p = — 1 to p' = 1, = 2vp J — |(r+r') + (r— r')j dr, 

— when the attracted point is without, and + when it is 
within the shell, 

= !Z r 2( a r '* dr ' = ^ } 
r J 0 3r 

when the point is without the sphere. 

When the point is within the sphere, the part of V for the 
shells which enclose the point 

= 27 rp J 2 r dr — 27 rp (a a — r*) : 
and the part of V for the other shells of the sphere 

4777 ) r ,, , , 4 

= ~r~J T dr = B 7rpr ' 

Hence V = for an external particle, 


2 

V = 2irpa 7 — - irpr 1 for an internal particle. 
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Prop. To find the attraction of a homogeneous body, differ- 
ing little from a sphere inform, on a particle without it. 

42. Since the attracted particle is without the attracting 
mass, we must expand V in a descending series of powers of r, 
and shall therefore use the first of the expressions for V in 
Art. 40. Let the mean radius of the body = a ; and let 
a (1 + y) be the variable radius, y being a function of ff 
and <o , and its square being neglected. 

Then, for the excess of the attracting mass over the sphere 
of which the radius = a, effecting the integration with respect 
to r' from r = a to r — a (1 + y), the value of V 

-pfX {; + + - +£*.+•••} 

But if y, the same function of p and o> that y is of f and 
a', be expanded in a series of Laplace’s Functions, viz. 

r 0 +F 1 +... + r i +..., 

then the theorems of Art. 26 and 32 show that 


2 i+ 


47 r 


- f 1 fy'Pffp'do' = P r Y/Pffp'da,' = F ( . 

J Q J __yj 0 


Hence the value of V for the excess over the sphere becomes 

— 47r P a . 3 ( y +— Y + + - Y+ 1- 

and the part of V for the sphere, rad. = a, is 

y_ iTT P &* 

3 r ' 

Hence for the whole mass 


3r 


(2 i+ 1) 




This is the first example in which we see the great value of 
the properties of Laplace’s Functions; they here give us at 
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once the integrals involved in our expression for V, in terms 
of the equation to the surface of the attracting mass, without 
integration. 

From the expression for V the attraction can be immediately 
found by the formula of Art. 20. Thus 

dV 

attraction = — 5— 
dr 


47rpa 3 i 47rpa 3 


\Y + - Y + 
1 0 3r 1 


(f + 1 ) a y 1 
(2t’+ l)r‘ 1 + 


Prop. To find the attraction of a homogeneous body, differ- 
ing but little from a sphere , on a particle within its mass. 


43. We must in this case expand V in an ascending series 
of powers of r; and shall therefore take the second of the 
series of Art. 40. By proceeding as in the last Proposition, 
we find that the part of V which appertains to the excess over 
the sphere 

= p [ f {a 2 P 0 + &rP x + ... + -ptP t + ...} y'd/i dm', 

J-iU a 




Adding to this the part of V which appertains to the sphere of 
2 

radius a, viz. 2irps? — - 7 rpr*, for the whole mass, 

O 

F=2 7 rpa 2 -|7rpr 8 +47rpa ! {F 0 + ^ Y t + ...+ Y l+ „.}. 

dV 

And the attraction = — 7— 

dr 



Y i + . 


ir 


(2»+l) a 1 





We can show that by properly choosing the value of (a) and 
the origin of the radius of the surface we can make Y 0 and Y 1 
disappear from the above formulae. 
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PROP. To show that hy choosing a equal to the radius of the 
sphere of which the mass equals that of the attracting body we 
cause l a to vanish, and by taking the centre of gravity of the 
body as the origin of the radius vector, we cause Y x to vanish. 


44. The mass of the body 

= p ff f^drd^-p f J* i*dpda>, 

where r is the radius vector of the surface of the body, and 
= a (1 +y) suppose. Putting this for r mass of body 

= mass of sphere (rad. = a) 4- pa 8 f f ydpdco 

J -i*' 0 

= mass of sphere + pa 8 f f F 0 , by Art. 26, 

^ -1*' 0 

= mass of sphere + 4 77-pa 5 F 0 . 

If then a be taken equal to the radius of the sphere of which 
the mass equals the mass of the body, F 0 = 0, as was stated. 

Again, let xy z be the co-ordinates to the centre of gravity 
of the body, M its mass : the co-ordinates to the element of 
which the mass is — pr 2 drdpdo> are 

rVrr p 1 cos to, rVf — p l sin < 0 , and rp ; 

tf r 


.•. M.x = J J j pr 3 Vl — p, a cos a>drdpdca 

1 r 1 r®*- 

= - 1 I pr 4 \/l — p? cos codpda, 

_i J 0 

M. y = f f f pr 8 vT — fi 2 sin co dr dfidco 

J oJ -I J 0 

2 r 1 r 2r 

= - I pr 4 Vl — p* sin codpdw, 

* J -1 J 0 
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M . z= J J J pr 3 pdrdpdco = ^ J J r *pdpda>; 

putting r = a (1 + y) = a (1 + Y 0 + Y, + ... + Y t + •••)> an( l °b- 
serving that VT cos od, Vl — p? sin co, and p satisfy 
Laplace’s Equation, and are of the first order, we have by 
Art. 26, 

M.x = p&*l j Y t Vl — p* cos codpdco, 
r i r2ir 

M.y = pa* I Y t Vl — pd sin co dpdco, 

J -l j 0 

M.z=p&*( f Y l pdpdo). 

But Y l , being a function of p, Vl — pi 1 cos co, and Vl — /a* sin co 
of the first order, is of the form 

A V 1 — p 2 cos c» + -5 ^ 1 — /a 2 sin cd + Cp ; 

_ 4 _ 4 4 

M. x = - 7rpa.*A, M.y =-~irpa,*B, M. z = - n-pa,* G. 

o o o 

Hence if we take the origin of co-ordinates at the centre of 
gravity and therefore x = 0, y = 0, z = 0, we have A = 0,B=0, 
C = 0, and therefore Y t = 0, as stated in the enunciation. 

Prop. To find the attraction of a heterogeneous body upon 
a particle without it; the body consisting of thin strata nearly 
spherical, homogeneous in themselves, but differing one from 
another in density. 

45. Let a' ( 1 + y') be the radius of the external surface of 
any stratum, a' being chosen so that 

y' = Y'+ Y; + ...+ F/ + ... (Art. 44). 

Since the strata are supposed not to be similar to one 
another, y is a function of a' as well as of p' and co'. Let 
p be the density of the stratum of which the mean radius 
is a. Now the value of V for this stratum equals the differ- 
ence between the values of V for two homogeneous bodies of 
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the density p and mean radii a and a' — da'. But for the 
body of which the mean radius is a' (Art. 42) 

4Va' 3 4 vpa' (a' , «" , \ 

3r r |3r I + # " + (2* + l)r‘ 

Hence for the stratum of which the external mean radius 
is a, 

and therefore for the whole body, 

v - t/V f' + v (£ r.'+ ••• + *+••■)} "• 

From which the attraction is easily deduced. 

Prop. To find the attraction of the same body on an in- 
ternal particle. 

46. Let r = a{l+y) be the radius of the stratum in which 
the attracted particle lies. Then for the strata within the 
surface of which the radius is a (1 +y), we have 

f= £]/{*+£(£ r ; + ... + 17+-.)}*. 

But for a stratum external to the particle we have by 
Art. 43, 

V = 17 + ••• + p *7+-}*. 

Consequently for the whole body, 

-r />' {*' + B (£ r -' + - + prnp Y - + •••)} *■ 

+ 4-/Vf'*'+^(T r -'+- + F ‘' + -)K 

From this the attraction is readily obtained by differ- 
entiating with respect to r. 








CHAPTER IV. 


ATTRACTION OF BODIES NEITHER SPHERICAL NOR SPHEROIDAL, 

NOR NEARLY SO. 

47. The methods which have hitherto been given enable 
us to find the attraction of the Earth and other bodies of our 
system considered as a whole. But, taking the Earth as our 
example, the surface is irregular, and neither exactly spherical 
nor spheroidal. We ought, therefore, to be able to calculate the 
effect of these irregularities, and with this view the present 
Chapter is added to what has gone before. High Table-lands 
may very materially affect the position of the plumb-line in 
some places. Enormous irregular mountain masses, like the 
Himmalayas, may do the same. Their effect ought, therefore, 
to be carefully estimated, as all instruments which are fixed 
by the plumb-line or spirit-level must be affected by such 
irregularities. 

^ Prop. To find the attraction of a slender prism of matter 
on a point in the line drawn to one of its extremities. 

48. Let AB be the prism, C the attracted point, P any 
element of the prism, AP—r, 

M the mass and l the length 
of the prism, AC=a, BC= h, 

PC = y, angle PA C = 9. 

Then the mass of the ele- 
ment at P= M y . 

Attraction of element at P 

n nr & 1 

on C=M- r -j. 

1 f 

Ditto in direction 

CA = M% cos PC A 

i y 
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_ Mdr a — r cos 0 

~~i 7~ 


if = a a + r* — 2 ar cos 0, 


Mdr a sin* 0 — cos 0 (r — a cos 0) 
l ’ (a a + r a — 2 ar cos 0 }* * 


attraction of whole prism 


M r — a cos 0 + a cos 0 
l a fa 1 + r a — 2 ar cos 0 


from r = 0 to r = l, 


JM l M 

cd Va a + Z* — 2 al cos 0 

As this is symmetrical with respect to a and b, it shows 
that the particle is attracted equally towards the two extre- 
mities of the prism ; and that therefore the resultant attraction 
acts in a line bisecting the angle which the prism subtends at 
the attracted point. 

Prop. To find the attraction of a slender pyramid of any 
form upon a particle at its vertex ; and also of a frustum of the 
pyramid. 

49. Let l be the length of the pyramid, a the area of a 
transverse section at distance unity from the vertex; r the 
distance of any section ; ar a is its area ; p the density of the 
matter : then ar^pdr is the mass of an element of the pyra- 
mid, and this divided by r a is its attraction ; 

.•. attraction of pyramid on vertex = apdr = apl. 

Jo 

If d is the length of any frustum of the pyramid, and l—l'+d, 
then 

attraction of pyramid, length T, = apT ; 

attraction of frustum = apd. 

It is observable that this is quite independent of the distance 
of the frustum from the vertex ; and therefore all portions of 


< 
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the pyramid of equal length, any where selected, attract the 
vertex equally. 

50. C()R. Suppose the angular width of the pyramid to 
be /3 and to remain constant, while the angular depth varies ; 
and let k be the linear depth of the transverse section of the 
base ; then afilk is the area of the base ; and the attraction of 
the whole pyramid on the vertex = p/3k. Hence, all pyramids 
having the same angular width and the same linear depth at 
the base attract their vertex alike, whatever their lengths be. 

Prop. To find the attraction of an extensive circular plain 
of given depth or thickness upon a station above its middle 
point. 


51. Let t be the thickness or depth ; h the height of the 
particle from the nearer surface, c the radius, r the radius of 
any intermediate elementary annulus of the attracting mass, 
z its depth. The several elements of this annulus of matter 
will attract the particle towards the plain equally. Hence 
attraction of the particle 




const. — 


Vr 2 +(/<+z) s 


dr 


: - r l 

dr 2 + ii* vV* + (h + <)*) 


= 2 irp {V c 2 + A 2 — h — Vc 2 + (A + t) s + h + £} 


= 27 rp i 

5-V7+A ? j(i + 

no 

II 

i h + \t ^ 

Vc 2 + A 2 


2ht + A* 
c 2 +W 


52. If the plain be of infinite extent, the attraction equals 
27 xpt ; and this remarkable result is true, that it is independent 
of the distance from the plain. The same will be the case if 
the height of the station above the middle of the attracting 
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mass below, that is, h + ^t, be so small that it may be neg- 
lected in comparison with the distance of the station from the 
furthest limit of the plain. 

53. Ex. Suppose the height of the station above the 
middle of the mass below, that is, h + \t, is £ a mile and c 
10 miles. Then the second term within the brackets is less 
than 0'05, and the attraction is very much the same as if the 
plain were unlimited in extent. 

54. Cor. The result of this Proposition when the plain 
is unlimited in extent might have been foreseen from the 
result in the previous Proposition regarding the attraction of 
the frustum of a pyramid. Conceive an infinite number of 
slender pyramids to be drawn from the station intersecting 
the attracting plain ; they will cut out of it an equal number 
of frustra, and the cosines of the angles they make with the 
perpendicular to the plain will be the thickness divided by 
the lengths of the frustra. But the attractions of the frustra 
are proportional to their lengths, and independent of the 
distance from the attracted point : (see Art. 49). Hence the 
resultant attraction of the whole will depend solely upon the 
thickness or depth of matter constituting the plain. 

PROP. To find the attraction of a rectangular mass , of 
small elevation compared with its length and breadth, upon a 
point lying in the plane of one of its larger sides. 

55. Let the attracted point be the origin of co-ordinates ; 
the axes of x and y parallel to the long edges of the tabular 
mass, the axis of z being measured upwards. Let x’y'z be 
the co-ordinates to any point of the mass : xy co-ordinates to 
the nearest angle, XY to the furthest angle, H the height of 
the mass ; p the density, supposed the same throughout. 

Then pdx'dy'dz is the mass of the element; and the 
height being small, we may suppose the element projected on 
the plane of xy. Hence the whole attraction parallel to x 

m " x'dxdy'dz' _ [ r x'dx'dy 

, i^ + yT~ p JJ sl^+yl 1 
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To simplify the formula put 

L- 


— = tan 0,, = tan 6 i , ~ = tan 0 3 , - = tan 

•C A A 3? 



F _ 1+ sin 0, 
a; cos 0, 


= tan (45° + 


and so of the rest. Hence, since 0'434 is the modulus of 
common logarithms, 


attraction = 


Pg 

0-434 


log tan (45° + ^ 0J + log tan (45° + ^P a ) 

- log tan (45° + £ d^j - log tan (45° + £0jj-, 


which gives a remarkably simple rule for finding the attraction 
parallel to x : that parallel to y can be found in like manner. 


It is easy to show, that if the density be half the mean 
density of the earth, that is, about the same as granite, g be 
gravity, the radius of the earth = 20923713 feet, and II be 

expressed in feet, the coefficient above = • 


This equals gH tan • Hence, since the tangent of 

deflexion of the plumb-line caused by the attraction equals, 
by the parallelogram of forces, the ratio of the attraction to 
gravity, and the angle is very small, 

Deflexion of plumb-line caused by the Tabular Mass 
parallel to the axis of x 
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= 7^9 H | lo S ten (45° + + log tan (45° + §<?„) 

- log tan (45° + £0 S ) - log tan (45° + £0 4 )j . 

It is evident that the Tabular Mass may be partly below 
and partly above the plane of xy , so long as the height or 
depth is not so great that its square may not be neglected in 
comparison with the square of the distance from the attracted 
point. In this case H is the sum of the height and depth, 
above and below the plane of xy. 

56. Ex. 1 . The co-ordinates to the nearest and furthest 
angles of a tabular block of rock measured from the attracted 
point are 3 and — 16, 40 and 30 miles, and the height of the 
mass from bottom to top is 628 feet. Show that the deflexion 
oT the plumb-line parallel to the shorter side of the parallelo- 
gram = 3''T72. 

Ex. 2 . A table-land 1610 feet high, commencing at a 
distance of 20 miles from Takal K’hera, near the Great Arc 
of Meridian in India, runs 80 miles north, and 60 miles to the 
east and 60 to the west. Find the deviation of the plumb- 
line at that station. It is about 5" ; so considerable as ma- 
terially to affect the Survey operations, and to have rendered 
it necessary to abandon that place as a principal station. 

In cases where the attracting mass is near, it is necessary to 
cut it up into prisms and calculate the effect of each separately 
and add the results. Examples of this are seen in the cele- 
brated case of Schehallien, and more recently in the calculation 
of the deflexion at Arthur’s Seat, Edinburgh, by Lieut. 
Colonel James, Superintendent of the Ordnance Survey. See 
Philosophical Transactions for 1856, p. 591. 

57. The irregular character of the surface of the Earth 
over large tracts of country, consisting of mountain and valley 
and ocean, may in some instances have a sensible effect, by 
presenting an excess or deficiency of attracting matter, upon 
the position of the plumb-line, in such a way as to derange 
delicate Survey operations. Hindostan affords a remarkable 
example of this, as the most extensive and the highest 
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mountain-ground in the world lies to the north of that con- 
tinent, and an unbroken expanse of ocean stretches south down 
to the south pole. Both these causes, by opposite effects, 
make the plumb-line hang somewhat northerly of the true 
vertical. 

In the following Propositions a method is laid down for 
calculating the attraction of an irregular superficial stratum 
of the Earth’s surface, and making it depend altogether upon 
the contour of the surface. The method pursued is this : A 
law of geometric dissection of the surface is discovered which 
divides it into a number of four-sided spaces, such that if the 
height of the attracting mass were the same in them all, they 
would all attract the given station exactly to the same amount, 
whether far or near. In this case it would be necessary only 
to calculate for one space, then count the number of spaces in 
the country under consideration, and the final result is easily 
attained. The country being supposed irregular, the heights 
in the spaces will not be all alike. The principle, therefore, 
should be stated thus, that the attractions of the masses on 
the several compartments are in proportion to their mean 
heights. These mean heights are known by knowing the 
contour of the country. 

-Prop. To discover a Law of Dissection of the surface of the 
earth into compartments, so that the attractions of the masses of 
matter standing on them, upon a given station, shall he exactly 
proportional to the mean heights of the masses, he they far 
or near. 

58. Suppose a number of great circles to be drawn from the 
station in question to the antipodes, making any angle $, eacli 
with the next, thus dividing the earth’s surface (which we 
may in this calculation suppose to be a sphere, without incur- 
ring any sensible error) into a number of Lunes. Then, with 
the station as centre, describe on the surface a number of cir- 
cles, at distances the law of which it is our object now to 
determine, dividing the whole into a number of four-sided 
compartments. 

We will begin by calculating the attraction of a mass of 
matter, standing on one of these compartments, at a uniform 
height throughout, upon the station in a horizontal direction, 
p. a. 4 
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Let a and a + <f> be the angular distances from the station of 
the two circles bounding this compartment ; h the height of 
the mass ; 0 the angular distance along the surface of an 
elementary vertical prism of the mass ; a the radius of the 
earth ; i/r the angle which the plane of 0 makes with the plane 
of the mid-line of the lune, and in which latter plane the re- 
sultant attraction evidently acts. The area of the base of the 
prism = a" sin 0dyjrd0. 

Since the height of prism ( h ) is supposed very small, the 
distances of its two extremities from the station may be taken 
to be the same, and = 2 a sin \0. Its attraction along the chord 
of 0, by Art. 48, 


pcfh sin 0d0d\fr 
4a* sin* £ 0 


Attraction along the tangent to 0 = 


ph sin 0d0dyjr 
4 sin* ^ 0 


cos i 0 ; 


.-. attraction along the tangent to the mid-line of the lune 
_ ph cos* $0d0dyjr 


2 sin 1 0 

attraction of the whole mass 


cos yfr; 


phC^riP cos* 0 

2 J a J- J/3 8in£0 


d0 . cos yfrdifr 


=* ph sm Ap — — f a d0 

K 2 ] a sin \0 

= 2 ph sin ^/S {log, + cos £ (« + (f>) — cos | a} 


“ ** Sin if> f l0 8. + ™ (4»+ ») ' *> H) 


= 4 ph sin £/9sin l<f> 


{sin (£a + i<£) 


-sin(|a + l<p) 
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neglecting only the cube and higher powers of sin 

= ph sin . ^ C03 (^ a +4 <ft) ne gi ec ti n g & c . 
r sin (ia + t^>) b • 


The law of dissection we shall choose will simplify this ; 
for we are to assume such a relation between <f> and a that the 
expression in <fr may be constant, in order to make the attrac- 
tion the same for all compartments in which h is the same 
or varying as h where the heights of the masses standing on 
the compartments are different. As the value of the constant 
to which we equal the function of a and </> is quite arbitrary, 
we will assume it such that when a and (f> are small, $ shall 

— cl 4 

= A a In this case it = — 10 , = — . 

10 £« + ^« 21 


Hence 


<f> cos 8 (^g + 1$) _ 
sin (^x + l(f>) ~ 21 


defines the Law of Dissection. 


( 1 ), 


The attraction of the mass standing on the compartment, 
in consequence, 

4 

= — phBm%p-, 

an exceedingly simple expression. We may obtain it in 
terms of gravity as follows. Let p the density be the same 
as that of the mountain Schehallien, viz. 2.75; the mean 
density, according to Mr Baily’s repetition of the Cavendish 
experiment, being 5.66 ; g gravity ; a = 4000 miles. 


-nt 47t , .. 4?r 566 

Now g = — ax mean density = — a — — p ; 

6 o * i o 

. attraction of mass on any compartment 

"DsS^ =0 - 00000552MBini ^’ 


being expressed in parts of a mile. 


4—2 
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Since 0.000005523 = tan ( 1 ". 1392 ) ; 

.*. deflection of the plumb-line caused by this attraction 

= 1 ". 1392 h sin ( 2 ). 

The method of using this theorem is as follows. When the 
numerical values of the successive pairs of a and <f> are deter- 
mined by the solution of equation (1) giving the law of 
dissection, lay them and the lunes down on a map of the 
country the attraction of which is to be found. It will thus 
be covered with compartments. After examining the map, 
write down the average heights of the masses standing on all 
the several compartments of any one lune ; add them together, 
multiply the sum by 1".1392 sin \ft, and equation ( 2 ) shows 
that we have the deflection caused by the mass on the whole 
lune in the vertical plane of its middle line. Multiply by the 
cosine and then the sine of the azimuth of that middle line, 
and we have the deflections in the meridian and the prime- 
vertical. The same being done for all the lunes, and the 
results added, we have the effects in meridian and prime- 
vertical produced by the whole country under consideration. 


59. Cor. That the mass on each compartment will attract 
as if collected at the middle of the mid-line appears from what 
follows. Let © be the distance at which the matter may 
be concentrated so as to produce the same effect as the actual 
mass. Now the area of the compartment 


favt> rip r°*t> 

= I a* sin ddddty — a? ft sin Odd 

J a J J a 


= a* ft {cos a — cos (a + <£)} = 2 a' ft sin %<f> sin (a + 


Then by the last Proposition, 

ph sin \ft = attraction, by hypothesis 

= 2a-g»in^,m(a + W 
r 4a* sm £© 2 ’ 


Digitized by GoogI 



DIMENSIONS OF THE COMPARTMENTS. 


53 


sin 8 _ -|y9 sin ^(f> sin 2 (^g + j<f>) 

cos^@~ sin^/S %<p cos (^a + ' 

Hence, if /9 be not taken larger than 30°, and since <f> is 
always small, this gives 0 = a + which coincides with the 
middle of the mid-line of the compartment. 

Prop. To calculate the dimensions of the successive compart- 
ments from the law of dissect ion. 

60. For this purpose we should solve the equation of last 
Proposition, viz. 

<f> co s 8 ($a + j<f>) _ £ m 

sin(£a + £<£) 21 ' 

But this cannot be done. We must therefore approximate, 
which will equally well suit our purpose. In order to afford 
a test of the values we arrive at the equation may be written 
under the following form. Putting the angle <f> for the 
arc (f>, 

£ 180 sin ($a + j<f >) . 

™ 21 it cos 2 (|a + £(/>) * 


(11.0379639 

or = log'V + log sin (|a + 

l- 2 log cos (£a 4- l<j>) 


(3). 


Equation (1) can be solved by expansion so long as a and 
<p are not too large. 

It gives 

= lr( a+ l){ i+ i(l + l)} ; 
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= ^(1 + 0.2411a*), 
or, if a be expressed in degrees, 

“ A ^ + °- 000073a ^ * (4). 


Let a^jOg ... ... be the successive values of a and <f> 

for the several compartments of a lune. These are connected 
by the following relations : 

«* =«! + &> «a =za » + < P»> 

Suppose that for the first compartment a 1 = 0°.75* ) then 


* This particular value is here used became the calculations following are 
taken from a Paper, by the author, in the Philosophical Transactions for 1855, 
upon Himmalavan Attraction, and three-fourths of a degree is about the dis- 
tance of the nearest hills from the northern station of the Great Indian Arc of 
Meridian. Any other value for a, might have been taken. The results above de- 
duced are perfectly general, and are applicable to any other similar problem. If 
compartments with elevations or depressions occur in the map, nearer to the station 
* than three-fourths of a degree, or about 52 miles, their dimensions can easily be 
calculated backwards, i.e. towards the station, by the help of the formula 

11. 

<P-- = a --ti i since = 

The succeeding values, reckoning from a, and 0 t are here given, as they may be of 
use for reference. 


a, = 0°.75 

a 0 0 .68 
a_, = 0 .62 

a-j = 0 .56 
a-j = 0 .51 
a_4 = 0 .46 
«_ 5 = 0 .42 
= 0 .38 
a — , = 0 .34 
0-8 — 0 ,31 


0, = 0°.07S 

0o =0.068 
0-, = 0 .002 
0-3 = 0 .056 
0-3 = 0 .051 
0-4 = 0 .046 
0- s = 0 .042 
0_s = 0 .038 
0-7 = 0 .034 
<p-e = 0 .031 


o -8 = (>‘>.28 

u_iq = 0 .2o 

a-,, = 0 .23 
a- ,8 = 0 .21 
a-jj = 0 .19 
n— 14 ” 0 .17 

0-15 = 6 .15 

a-,8 = 0 .14 
a_, 7 = 0 .13 
he. 


0-s = 0°.028 

0— io = 0 .025 
0— ti ~ 0 .023 

0-u = 0 .021 

0—13 = O .019 
0—14 = 0 .017 
0—15 = 0 .015 
0—is = 0 .014 
0— ,7 = 0 .013 
he. 
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by (4) 0, = 0°. 075 ; therefore = 0°. 825 : and by proceeding in 
this way we obtain the following pairs of values by this 
formula. 


a, =0°.75 0, =0°.075 
c s = 0 .83 fa = 0 .083 

0 3 =0.91 0, =0.091 

0 4 = 1 .00 fa = 0.100 
o , = 1.10 < j b s = 0.110 
a 0 = 1.21 fa = 0.121 
a, =1.33 <£, =0.133 
o 8 = 1 .46 <£„ =0.146 
a 9 =1.61 fa =0.161 
°io = 1 *77 fa 0 = 0.177 
a u =1.95 0 U = 0.195 
a 12 = 2 .14 <k 2 =0.214 
a ls = 2.35 <£,,= 0.235 
a 14 = 2 .59 <£„ = 0 .259 
a, 5 = 2 .85 • <£„ = 0 .285 
o,„ = 3.13 0„, = 0.313 
a„ = 3 .45 <£,, = 0.345 
c 18 = 3.79 <£,, = 0.379 
o, 9 = 4.17 <£, 9 = 0.417 
a M = 4 .59 fa m = 0.459 
a si = 5 .05 <£„, = 0 .505 


o^ = 5° .55 

033 = 0°.555 

033=6 .11 

033=0.611 

o M = 6 .72 

0^ = 0 .672 

033 = 7 .39 

0.; = 0.739 

033=8 .13 

0m = O >813 

Oj, = 8 .94 

03, = 0.894 

o a = 9 .83 

0 2a = O.983 

0^= 10.82 

0^=1.089 

«»= H-91 

030=1 -202 

o 31 = 13.11 

0 31 =i.326 

033=14.43 

033 = 1-462 

0^= 15.99 

0 33 = 1 .62O 

03,= 17.61 

0, = 1.800 

033=19.41 

<5^=1 .992 

033 = 21.40 

0 m = 2.211 

o 87 = 23.61 

03, = 2.456 

o M = 26.06 

0,3 = 2 .732 

o w = 28.79 

0 W = 3.O54 

c 40 = 31.84 

040 = 3 .419 

o 41 = 35.26 

0 4J = 3 .600 

o 4a = 38.86 

0,3 = 4.314 


That these results derived from the approximate formula 
(4) are thus far correct, we gather from the fact that the last 
pair, viz. a =38°.86 and <£ 42 = 4°.314 sufficiently satisfy the 
test (3) when substituted. Beyond this pair, we cannot 
use (4), but must solve equation (1), or rather (3), by trial. 
This leads to the following pairs of values stretching to 
the antipodes. 


0,3 = 43°. 17 0 4S = 4°.98O 
o 44 = 48 .15 0„ = 5.783 
0^ = 53.93 <£ 45 = 6.800 
o*. = 60 .73 <£*, = 8.210 


o 4 ,= 68°.94 0„ = 10°. 330 
o w = 79.27 <£*,= 14.030 
o 49 = 93.30 0,3 = 23.380 
o w = 116 .68 ^ imperfect. 
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61. The formulas here deduced may be applied to find the 
effect on the plumb-line of any mountain-region, or hollow (as 
in the case of the ocean), so long as the angle subtended at 
the station by any part of it is such as to allow its square to 
be neglected 

In the Philosophical Transactions for 1855 and 1858-9, the 
author has applied these principles to find the effect of the 
Himmalayas and the mountain-region beyond them on the 
plumb-line in India, and lias found that the meridian deflection 
caused in the northern station of the Great Arc of Meridian 
(lat. 29° 30' 48", and long. 77° 42') is nearly 28", as far as the 
data regarding the contour of the mass can be ascertained; 
and that the astronomical amplitudes between that and the 
next principal station (lat. 24° 7' 11"), and between that and 
the third (lat. 18° 3' 15"), are diminished by the quantities 15". 9 
and 5". 3. He has also shown that the meridian deflection be- 
tween the first and third of these stations varies very nearly 
inversely as the distance from a point in the meridian in lati- 
tude 33 s 30'. 

62. The effect of the deficiency of matter in the Ocean 
south of Hindostan down to the south pole is also calculated, 
upon an assumed but not improbable law of the depth, and 
found to produce a meridian deflection northwards at the three 
stations specified of about 6", 9", 10". 5 respectively; and 
19".7 at Cape Comorin. 

63. It is possible that, the superabundant matter in 
mountain-regions having been heaved up from below, there 
may be a deficiency of matter below the mountains which 
would under certain circumstances have the tendency of 
counteracting their effect on the plumb-line. This Mr Airy 
has suggested in a Paper in the Philosophical Transactions 
of 1855, on the hypothesis that the deficiency is immediately 
below the mountains close to their mass. Upon the supposi- 
tion that the mountains may have drawn their mass from 
the regions below through a considerable depth, by an exten- 
sive and small expansion of the matter in those lower regions, 
the author has calculated the modifying effect on the plumb- 
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line in the Transactions for 1858-9. This has brought to light 
the fact, that a trifling deviation in the density from that 
required for fluid-equilibrium, if it prevail through extensive 
tracts, may have a sensible effect upon the plumb-line. The 
following Proposition, with which we shall close this Chapter, 
will show this. These questions, in themselves interesting as 
problems in Attraction, become still more so, as we shall see, 
in the determination of the Figure of the Earth. 

Prop. To find the effect on the plumb-line of a slight but 
wide-spread deviation in density in the interior of the earth, 
either in excess or defect, from that required by the laws of 
fluid-equilibrium. 


64. Suppose vertical lines drawn down through the four 
angles of any compartment to a depth d, and a surface uniting 
the four extremities drawn, so as to form the frustum of a 
pyramid of which the vertex is in the centre of the earth : 
draw also a vertical line of length d through the mid-point of 
the compartment. Suppose the height of the matter standing 
on the compartment to be uniform and equal to one mile. 
Let the several vertical prisms of which it consists be con- 
ceived to be distributed downwards uniformly through the 
depth d, the density of this lengthened prism will be less 
than that of the superficial rock in the ratio of 1 : d. Let u 
and v be the distances of the extremities of this long prism 
from the station. Then the attraction of the short prism 


along the chord of the surface 


mass 
u * 


, and (by Art. 48) that 


of the longer = 


mass 

uv 


Hence along the horizontal line at the 


station 


attraction of slender prism 


attraction of the prism at surface 


u 

v 


Now in Art. 59, it has been shown that the attracting mass 
on any compartment may be considered concentrated in the 
mid-point. Much more may this be done with the horizontal 
layers of the frustum which are not of larger dimensions than 
the compartment, and are farther off from the station. Hence 
if anti v, be the distances from the station of the extremities 
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of the vertical line d through the middle point of the mid-line 
of the compartment, the attraction of the mass on the com- 
partment, and the deflection caused by it, must both be 
diminished in the ratio of u. to v t to find the effect of the 
same mass distributed through a depth d. Suppose the masses 
(one mile high) on n compartments of any lune are thus 
distributed ; then by Art. 58, formula (2), 


Deflexion = l". 1392 sin 


* + * + 



If £ = 30°, the coefficient = 1". 1392 x 0.258 = 0". 294. 


65. We will take an example. Let the width of the lune 
/3 = 30°, and let the 21 compartments from a, to (see Table 
in Art. 60) be included. This will be a tract of country 
5°. 55 — 0°.75 = 4°. 8, or 334 miles in length, and the breadth 
at the mid-point will = sin ^ (5°. 55 + 0°. 75) x the length 
of 30°= 0.055 x 30 x 69.5 = 114 miles; and, by spherical tri- 
gonometry, the area is, in round numbers, 38,500 square miles. 
We will take three examples of depth which (for convenience 
of calculation) we will express in the length of degrees, viz. 
3®, 6°, and 9°; which nearly equal 208, 417, and 625 miles. 
The vertical thicknesses of these three divisions of the frustum 
are each = 3° = 208 miles. The widths, however, parallel to 
the horizon grow less in passing downwards. But owing to 
the convergency of the radii bounding the elementary prisms, 
the density increases in the distribution of the matter in 
exactly the same proportion that the area of the horizontal 
section diminishes. The amount of matter in the three divi- 
sions is therefore the same, and we may consider the volumes 
the same, and each equal to 38,500 x 208 = 8,008,000 cubic 
miles = 3-100, 000th parts of the volume of the whole earth. 

Now since the greatest of ... u n is less than 6°, we may- 
take the arc for the chord without sensible error. Then, with 

• ■ u 

respect to all these quantities, - = the cosine of the angle of 
which is the cotangent. This enables us without difficulty 
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with the help of a Table of cosines and co-tangents to form 
the sum of the series in the last Article. The values of u are 
the first 21 values of a in the Table in Art. 60. 


d= 

= 3° 

d-- 

= 6° 

d = 9° 

ej s 

II 

8 

p+ 

u 

cos — — 

V 

^ = cot 
a 

u 

cos = - 
V 

u 

= cot 
d 

u 

cos = - 

V 

0.250 

0.242 

0.125 

0.124 

0.083 

0.083 

0.277 

0.267 

0.138 

0.137 

0.092 

0.092 

0.303 

0.290 

0.152 

0.150 

0.101 

0.100 

0.333 

0.316 

0.166 

0.164 

0.111 

0.110 

0.367 

0.344 

0.184 

0.181 

0.122 

0.121 

0.403 

0.374 

0.201 

0.197 

0.134 

0.133 

0.443 

0.405 

0.222 

0.217 

0.148 

0.146 

0.487 

0.438 

0.243 

0.236 

0.162 

0.160 

0.537 

0.473 

0.269 

0.260 

0.179 

0.176 

0.590 

0.508 

0.295 

0.283 

0.197 

0.193 

0.650 

0.545 

0.325 

0.309 

0.217 

0.212 

0.713 

0.581 

0.356 

0.335 

0.238 

0.232 

0.783 

0.617 

0.392 

0.365 

0.261 

0.253 

0.863 

0.653 

0.431 

0.396 

0.288 

0.277 

0.950 

0.689 

0.475 

0.429 

0.317 

0.302 

1.043 

0.722 

0.522 

0.463 

0.348 

0.329 

1.150 

0.755 

0.575 

0.498 

0.383 

0.358 

1.263 

0.784 

0.631 

0.534 

0.421 

0.388 

1.390 

0.812 

0.695 

0.571 

0.463 

0.420 

1.530 

0.837 

0.765 

0.608 

0.510 

0.454 

1.683 

0.860 

0.842 

0.644 

0.561 

0.489 

sums = 

11.512 

= 

7.101 

= 

5.028 

multiply by 





0''.294 = 

=3''.385 

= 

CO 

CO 

o 

*CN 

== 

1".478 


From this Table we gather, that the Deflections caused at 
the station by the superficial mass one mile thick, when dis- 
tributed uniformly through the depths 208, 417, 625 miles, 
are 3”. 385, 2". 088, 1".478. The densities of the matter thus 
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diffused in these three cases are about -^th, y^th, 7 J^th 
of the density of the superficial rock. If we multiply the 
above deflections by 2, 4, 6, we have the deflections caused by 
matter, of l-100th the density of superficial rock, distributed 
over the three depths, etjual to 6".770, 8". 352, 8''. 868. Re- 
taining the first, subtracting the first from the second, and the 
second from the third, we have the three deflections caused by 
a mass of 208 miles vertical thickness (occupying 3-1 0Q, 000 th 
parts of the volume of the whole earth and of density 1-100 th 
part of the density of the surface) the centre of which is at 
depths 104, 312, 521 miles: they are 6". 770, l".582, 0".516. 
We may finally change the comparison between the density 
of this space, 3-1 00,000 ths of the earth’s volume, in its three 
situations, with the density of surface, to a comparison with 
the average density of the earth itself at the several depths 
at which the centre of the space lies. 


If D be the density of the surface, a the radius of the earth, 
the usually received law of density of the interior, determined 
from the fluid-theory, is 


Density at depth d = 


2 aD . 


vSin 

a — a 



When d= 100, 300, 500 miles, this gives the densities 
1.14D, 1.43D, 1.7 ID. 

Multiplying the last angles by the ratios of these densities 
to D, we have finally the Deflections — caused by an excess 
or defect of matter, prevailing through a space equal to 
3-100,000 th parts of the volume of the earth, and 1-100 th part 
of the earth’s density at the centre of the space— equal 7''. 7, 
2". 3, and 0".9, the depths of the centre of* the space being 
about 100, 300, 500 miles. 

The form of the space in its three positions is shown in the 
diagram ; viz. DG, FI, UK; and O, P, Q are their centres. 
The particular form arises from the manner in which the mass 
is dissected, so as to make the calculation feasible. The result 
serves to show the kind of effect which slight but extensive 
variations from the density of fluid-equilibrium in the hidden 
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regions below may have upon the plumb-line : and we shall 
find the use of this when we come to consider the Figure of 
the Earth : (see Art. 98). 


66. Had the width DB been equal to the middle width 
at a, so as to make the boundaries BG, DE parallel, the effect 

A B 



would have been very much greater. Moreover the defect or 
excess in density which we have taken, viz. 1-100 th, might 
have been chosen larger, and the deflections proportionally 
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increased. For there are many kinds of rock, as granite, 
which differ so in density in the different specimens that the 
difference between the extremes is greater even than 1-1 Oth 
of the mean. And if this difference exists at the surface, it 
does not seem to be improper to suppose that great variations 
may exist also below, from the effect of the cooling down and 
solidifying of the crust, even much greater than 1-100 th. 
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CHAPTER I. 

THE FIGURE OF THE EARTH CONSIDERED AS A FLUID MASS, 
AND THEREFORE CONSISTING OF STRATA 
NEARLY SPHERICAL. 

67. After it was known that the earth is of a globular 
form, Newton was the first who demonstrated that it is not 
a perfect sphere. From theoretical considerations and also 
from the discovery that a pendulum moves slower at the 
equator than in higher latitudes, he arrived at the conclusion 
that its form is that of an oblate spheroid. This subject we 
propose to consider fully in the present Chapter, on the 
hypothesis that the Earth was a fluid mass when it assumed 
its present general form. The calculation is one of great 
difficulty, aud would indeed be impracticable did we not know 
that the figure differs but little from a sphere. 

As a first approximation we shall inquire whether a ho- 
mogeneous fluid mass revolving about a fixed axis can be 
made to maintain a spheroidal form according to the laws of 
fluid pressure. 

Prop. A homogeneous mass of fluid in the form of a 
spheroid revolves with a uniform velocity about an axis : re- 
quired to determine whether the equilibrium of the surface left 
free is possible. 

68. Let a and c be the semi-axes of the spheroid referred 
to three axes of rectangular co-ordinates, c being that about 
which it revolves: also let c 8 = o s (1 — e a ). The forces which 
act upon the particle (xyz) are the centrifugal force and the 
attraction of the spheroid parallel to the axes: these latter 
are given in Art. 12, and are 
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{Vl — e‘ sin 1 e — e (1 — e*) j x, 
2 ^- {Vl — e‘ 8 in T ‘ e — e (1 — e*)} y, 


{e — Vl — e 2 sin -1 e\ z. 
e 1 ‘ 


Let these be represented by Ax, By, Cz. Let w be the 
angular velocity of the rotation, then w 2 V ar* + y* is the cen- 
trifugal force of the particle (xyz), and the resolved parts of it 
parallel to the axes of x, y, z are w*x, ufy, 0 . Hence X, Y, Z, 
the forces acting on {xyz) parallel to the axes, are 

X = — [A — w 1 ) x, Y= — {B—w v )y, Z= — Cz. 

These make Xdx + Ydy + Zdz a perfect differential, and 
therefore so far the equilibrium is possible. 

The equation of fluid equilibrium gives 


: dp = Xdx + Ydy + Zdz 

= — (A — w*) {xdx + ydy) — Czdz ; 



= constant — {A — v?) {a? + y 1 ) — Cz *. 


At the surface p — 0, and therefore 

^ q °- (as 2 + y 2 ) 4- z 2 = const. 


is the equation to the surface; and this is a spheroid, and 
therefore the equilibrium is possible, the form of the spheroid 
being properly assumed. The eccentricity is given by the 
condition 


1 — e* = -, = 


w 


c 


u> Vl - e 1 . 1 - e 2 2 . s 

or - — = = — sin e — 3 — * — h -3 (1 — e~\ l sm e ; 

2 irp e 8 e 2 e 8 v ' 
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v? 1 — e a 

**• o b 3 — g — 

2irp e 


(3 — 2e s ) Jl — e 2 


sin 1 e = 0. 


Now observation shows that ~ = the ratio of the centri- 

*o«/ 

fugal force at the equator to gravity at the equator. Hence 

1 _ w 2 a _ w 2 _ \ 

289 1 7 rpa — w 2 a ’ ' ‘ 2wp 435 " 


By expanding in powers of e and neglecting powers higher 
than the second, because we know that the earth is nearly 
spherical, we have 


... 1 e 3 1 . 3 e 5 

sm e-e + 5 - + 



If e be the ellipticity, then 


e 


a — c 
a 


= 1 - Vl -e* = 



1 

232 ' 


This result is so much greater than that obtained by other 
methods, as we shall see, that it decides against our consider- 
ing the earth’s mass to be homogeneous. Indeed it is h priori 
highly improbable that the mass should be homogeneous, 
since the pressure must increase in passing towards the centre 
and the matter be in consequence compressed. 

p. a. 5 
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69. Another value of e, nearly = 1, satisfies the equation. 
But this does not give the figure of any of the heavenly bodies, 
since none of them are very elliptical. 

Since there are two values of e which satisfy the equation, 
it might be supposed that the equilibrium of the mass under 
one of these forms would be unstable, and, upon any derange- 
ment taking place, the fluid would pass to the other as a 
stable form. But Laplace has shown (MSc. Cedes. Liv. ill. 
§21) that for a given primitive impulse there is but one form. 
In fact it is easily seen that for a given value of w, the angular 
velocity, the vis viva of two equal masses, so different in 
their form as to have e small and nearly equal unity, must be 
very different, and that therefore the mass cannot pass from 
one form to the other without a new impulse from without 
being given to its parts. 

70. The relation between w and e in Art. 68, shows that 

as w alters e alters, and vice versft. By putting — = 0, we 

find the greatest value of to which is consistent with equi- 
librium. This after some long numerical calculations gives 

17197 

e = — - — , and time of rotation = 0 - 1009 day. 

27197 J 

71. Before proceeding to calculate the ellipticity on the 
hypothesis of the earth’s mass being heterogeneous we will 
take the following extreme case. The density increases as 
we pass down towards the centre. Suppose that at the centre 
it is infinitely greater than elsewhere: that is, suppose the 
whole force resides in the centre. The case of nature must lie 
between this hypothesis and that of the earth’s being homo- 
geneous. 

Prop. To calculate the ellipticity of a mass of fluid revolving 
about a fixed axis and attracted by a force residing wholly in 
the centre of the fluid and varying inversely as the square of the 
distance. 

72. Let M be the mass of the fluid; the other quantities 
as before ; 
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.. Z= — ^- + wx, Y=-^f + w*y, Z=--^-. 

Then the equation Xdx + Ydy + Zdz = 0 becomes 
M 

(ax£c + ydy + zdz) — w 1 ( xdx + ydy) = 0 ; 

••• ^*-y<i(^+jfl = 0; 

— f — (a* + y 2 ) = constant = C. 


Then ™ = 


289 if 


— = — w a 
a 


4 ^ = 290; 
wa 


G 1 at s + y* 


r Va 2 + y 2 + js a M 580 a 8 

By reversing this, squaring, expanding, and neglecting the 
square of , this is seen to be the equation to a spheroid. 

When x=0 and y— 0, then z — c\ when z — 0, »*+y* = a* ; 

1 G_ 1 = £ _ J_ 1 c 580 , 

c if ’ a if 580 a ’ a 581 ’ 

1 

••• e = -~. 


This value of e is too small (as we might have expected), 
as — is too large, to agree with the form deduced in other 

2iO& 

ways. 

Prop. To find the equation of equilibrium of a heterogeneous 
mass of fluid consisting of strata each nearly spherical, and re- 
volving about a fixed axis passing through the centre of gravity 
with a uniform angular velocity. 

5—2 
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73. Let XYZ be the sums of the resolved parts of all the 
forces which act upon any particle (xyz) of the fluid, parallel 
to the axes of co-ordinates, p the density at that poiut, p the 
pressure. Then the equation of fluid equilibrium is 

— Xdx + Ydy + Zdz. 

At the surface, and also throughout any internal stratum 
of equal pressure and therefore of equal density, in passing 
from point to point dp = 0. 

Hence Xdx + Ydy + Zdz — 0 

is the differential equation to the exterior surface and to the 
surfaces of all the internal strata ; the particular value assigned 
to the constant after integration determining to which surface 
the integral belongs. 

The following property belongs to all these surfaces. If ds 
be the element of any curve drawn on the surface through 
(xyz), and It be the resultant of XYZ; then the equation may 
be written 


X dx Y dy Zdz_ 

R + R ds+R ds = °’ 

which shows that the resultant force is at right angles to any 
line in the surface, and therefore to the surface itself at the 
point (xyz). 

The equilibrium will be the same if we suppose the rotatory 
motion not to exist, but apply to each particle a force equal to 
the centrifugal force caused by the rotation. The forces then 
acting on the fluid will be the centrifugal force and the mutual 
attraction of the parts of the fluid. Let V be the potential 
Art. 18) for this mass, then 

_dV _dV _dV 
dx ’ dy ’ dz 

are the attractions parallel to the three axes tending towards 
the origin of co-ordinates. Let w be the angular velocity of 
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rotation about the axis of z, taken as the fixed axis; w t .x 
and w i .y will be the centrifugal force at the point (xyz ) . Then 
the differential equation to the surface and the strata becomes 



w 

or constant = V+ — (of +y r ). 

/t 

Let r be the distance of the point (xyz) from the origin, 
and 6 the angle r makes with the axis of z, and cos 9 = ji: 
then x* + y* = r 3 sin 3 6 = (1 — y?) r 2 . Also let m be the ratio 
of the centrifugal force at the equator to gravity at the 

equator (° r 2^9)’ ^ et a ' t ^ ie mean ra dius of the stratum 

through (xyz) ; a the radius of the equator ; then 

„ M m>V 
m = w a — 0= — rr , 
a 2 M 


r* t 4 

and M= 47 r pa n da ' = - 7r$(a) suppose, 

J 0 o 

the strata being considered spherical because of the smallness 
of the numerator in the value of m ; 


3 m? 2 a 3 


m = 


2 4<7r <f> (a) 

w — — m 

3 a 3 


47 Tip (a) ’ 
and the equation becomes 

constant = V+ ^ m (1 — y?) r 2 

O £1 

rr , 4tT <)>(&) , 27 T </> (a) /I A , 

= f + t w ir r+ T m z’ J (3-> l ) r - 


this arrangement being made, because the second and third 
terms as they now stand, are Laplace’s Functions of the 
order 0 and 2. (See Art. 39 , Ex. 1.) 

Now since the mass is supposed to be fluid and the external 
surface nearly spherical, it follows that as the heavier parts, 


■s 
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which are all free to move, will sink through the lighter, and 
lie in layers, these layers will also be nearly spherical, other- 
wise there will be a greater pressure on one part than on an- 
other, and the equilibrium of the layer or stratum will not 
exist. We may therefore assume as a consequence of the 
fluidity and the form of the surface that the strata also are 
nearly spherical. 


By Art. 46, 

we have 









jLt 

'a' 4 

Y' 4- 

4_ 

a' 1 ** 

Y' 4. M 

■da! 

V ' j/1 

il ~T 

da'\ 

,3r 


• • i 

{2i+l)j* 


+ 4w [ p - 
J a 

|a’ + 

-1 
da ' 

far 

vT 

y;+. 

• • + 

r‘ 

(2»+l)o' M 

r/+...) 

jda'. 

In this put 









r = a (1 + Y 1 - 

4“ ••• 

Y t 

+ ...) 

and 

j 0 p'a^da' 

d 0 

= 0(a), 



as before. Then substitute this value of V in the equation 
to the strata and equate terms of the order ». (See Art. 35.) 

The constant parts give 


[dp 4w <f> (a) 47 r 

Jf=T — 


ma 


0(a) 


and the terms of the order i give 


4 >( a ) 

3 a 


o* f , d ( y; \ , , 

2i+lJ a p da' la'‘- 2 ) da ~ 0 ’ 


except when i — 2, in which case the second side is 

m (a) /I _ A 
6 a 8 \3 W* 

By this equation T, is to be calculated, and then the form 
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of the stratum of which the mean radius is a is known by the 
formula 

r = a ( 1 + + l s +...+ 1 , + 

Prop. To 'prove that Y { = 0, excepting the case of i— 2. 

74. Since F, and p are functions of a, they may be ex- 
panded into ascending series of the form 

Y,= Wa‘+..., p = D + D'a* + 

v here D is the density at the centre of the earth, and is as 
well as W and D' independent of a : s, n ... must not be 
negative, otherwise F, and p would be infinite at the centre. 

Now when these and the corresponding series obtained by 
putting a for a, are substituted in the equation of the strata 
in the last Article, and the first side arranged in powers of a, 
the various coefficients ought to vanish ; excepting when * = 2, 
because then the second side is not zero. We shall therefore 
substitute these series, and search for values of W and s which 
satisfy the condition. 

r 8/r 

, <f> (a) — 3 pa* da = Da 3 + a“ +s + . .. 

J 0 W + o 

After two easy integrations the equation of the strata 
becomes 

No valuq of s will cause these terms to vanish. The only 
apparent case is when t = 1, for then by putting 8 = * — 2 the 
part in the brackets vanishes : but in this particular case 
s = — 1, and is negative and therefore inadmissible. 

Hence the only way of satisfying the condition is by putting 
W= 0 ; this shows that F, has no first term, that is, that it 
has no term at all and is therefore zero. 

Prop. To find the value of F s , and to prove that the strata 
are all spheroidal , concentric and with a common axis. 
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75. The equation for calculating Y t is, by Art. 73, 



m a*6 (a) t\ 
= 6 


Suppose T s (and similarly Yf) is expanded in a series of 
powers of J — /i“ with indeterminate coefficients to be ascertain- 
ed by the condition, that they shall satisfy the above equation. 
These coefficients will be functions of a only, as it is clear from 
the right-hand side of the equation that &> does not enter into 
the value of Y t ; and T s consists of only one term, that in- 
volving the simple power of J Let it be e (J — /x s ), e 

being a small quantity of the order of in. Hence 

r = a (1 + £ (J — /**)}, fi = sin (latitude) = sin l 
= a (1 — § g) (1 + s cos* l), since £ is small. 

This is the equation to a spheroid from the centre, e being 
the ellipticity. The axis-minor coincides with the axis of 
revolution of the whole mass. Hence the strata are concen- 
tric spheroids, the minor-axes of which coincide with the axis 
of revolution of the whole mass. 

Prop. To obtain an approximate law of the density of the 
strata. 


76. By Art. 73 we have the following equation for calcu- 
lating the pressure on the stratum of which the radius is a, 
neglecting the small term, 

» 



p da 3 a 


+ 47 r 


I p'a'da'. 

J a 


Laplace has integrated this equation on the supposition 
that the change in pressure in descending through the strata 
varies as the change in the square of the density ( Mdmoires 
de V Institut, Tom. ill. p. 496). This law of compression 
differs from that of elastic fluids, in which the change in 
pressure varies as the change in density. The law used by 
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Laplace is h priori more probably true than the law of 
compression of elastic fluids, for the greater the density of 
tenacious and semi-fluid masses the greater must be the increase 
of pressure to produce a given increase of density. See also 
some remarks on this subject by Professor Challis in the 
Phil. Mag. Vol. xxxvill. The approximate truth of this law 
is, however, shown by the accuracy of the results to which it 
leads us. 


Putting, then, dp = \kd.p*, k being a constant, 

[ da' = k(p + constant) ; 

J 0 p da 

ka{p + const.) = 4 v / pa* da' + 47ra I p'a'da', 

J 0 •'a 

since <f> (a) = 3 p'a* da'. Differentiate* with respect to a ; 

*' o 

.'. ft +const.^ — iTrpa* + 47 t |" p a da! —K it pdJ—iir J p'a'da'; 

dJ .pa n _ • 4 7r . 

-jnr + f-pa=0, putting -r- = f ; 


pa = Q sin ( qa + B ) ; p = — sin ( qa + B). 

When a = 0, p = ® ^ B= 0, otherwise p would be 
infinite at the centre, which cannot be ; 

Q • 

.’. p = — sin qa, 

Q and q being unknown constants. 


• In order to explain how to differentiate a definite integral with respect to a 
quantity involved in the limits, let ff(x) = F (i) + const. ; 

.-. f b a f(x)dx=F(fi)~F(b)i 

d /*■ . , dF(a) ,, . d , dF(b) 

Ta sr =f{a)i Tb Jo fw* r =-~2r=-m- 
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Prop. To obtain an equation for calculating the ellipticity 
of the strata. 

77. Substitute e (J — p 1 ) for Y. and s' (J — fT) for Y t ' in the 
equation of the last Proposition but one, and we have, after 
dividing by § — p*, 

<f>(a) 1 [" , d . , 5 , , a* f a , de , , m a*(f>(a) 

3a 5 a 9 J 0 ^da K ' 5 ] a ^ da 6 a 3 

Divide both sides by a*, and differentiate with respect to a ; 
then multiply by a®, and differentiate again, and divide by the 

coefficient of -7-5 ; 

da 


jGpa 4 de f pa* 
(f>(a) da \ <f> (a) 


J*£ ( 6pa J de 


S=o. 


This may be put into another form. Multiply by <f> (a), then 

Cor. 1. By putting a = a in the first equation of this last 
Article, we have the following equation, which we shall find 
of use; 

/„ p i (“ V) *' - 1 a ’* <*> ( e - f ) • 

Prop. To find an expression for the ellipticity of the strata, 
with the law of density deduced in the last Proposition but one. 

78. In the equation of last Article put p = sin qa. 

Now <p(a) =3 J p'a*da =3 Q [ — cos qa + -3 sin qa}. 
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Als ° S = « {! C0S S» - J sin 3 a } - - ^ 4, (a), 
and our equation becomes 

jnH ~+ 2 *. <£(«) 6 = ~ S <#> («)•*• 

To integrate this put <f> (a) s = — a f a' f ax' da * ; 

a K Jo 

. d.<j>(a)e 2 f* , ^ , 1 f* , ,, , 

.. ”7 = i a I ax da +■- ax da: 

da a 3 J 0 J 4 aj 4 

. d\<f>(a)B 6 2 r , , 7 , 

da aV 0 J 0 <* J 0 « J 0 

x — ^J axda + ^ j a' j" ax' da? = 0. 

Multiply by a* and differentiate ; 

r a 

a* ^ + 2aa; — 3 aa? + j*a J a'x'da' = 0. 

Divide by a and differentiate, and then divide by a ; 
a?“a; „ 

55 + s*-°- 

The solution of this is 

a; + ( 7 j s sin (ja + _B) = 0, 

(7 and B being independent of a ; 

f a'x'da! = — Cqa cos (ja + 1 ?) + (7 sin (ja 4- 5 ) ; 

^ 0 

If 2(7 

(f> (a) e = ^5 j Ca* sin (gu + B) + -— a cos (ga + Z?) 

2(7 (7 (7 . 1 

— sin (ga + B) + — a cos (ga + B) — sin (ga + J 5 )| 

= C |(i - q * J) sin (ga + B) + -j^ cos (ga + -0)J . 
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In our case B = 0, otherwise the ellipticity at the centre 
would be infinite, as is easily seen by expanding g in powers 
of a. 


Hence, if we substitute for <f> (a) 


£ = 


3 Q tan qa — qa 

And e being the ellipticity of the surface, 


£ _ tan 17 a — qa. 
e tan qa — qa 




qa 




tan qa + — 
1 qa, 


This gives the law of decrease in the ellipticity of the strata 
in passing down from the surface to the centre. 

By Art. 77, Cor. 

»’«») («-§«•)- ■/ V TS («’■*')*>'= Q '/* («’■') id 

= <3 {a 4 e sin qa. + I a’V (sin qa — qa! cos ya') da) by parts. 

* 0 

Substituting for s' from the expression already found for e, 
integrating and reducing, the integral in this expression 

e (tan qa — qa!) sin qa, | 5 V — 15ya] 

q* (. 3 \ . 3 { 2 a tan qa, j 


S a 


30 


Also *(»)- 

Hence ~ = 
2e 


gV- «gV+ + (l - -? a -) ((SjV-lS - ^- 15 ? - a ) 
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Put 


— — = 1 — z to facilitate the calculation ; 
tan ga 


bin z (gV — 3 z) 

2 — gV + 3gVz — g*aV 



When this is calculated for the surface, we shall he able to 
find the ellipticity of any stratum we like by the ratio of £ to 
e already found above. 


Prop. To prove that the ellipticity of the strata decreases 
from the surface towards the centre. 

79. We assume that the density of the Earth increases 
from the surface to the centre. Let then p — D — Ed' + . . . , 
where E is positive : and £ = A + Ba m + .... Then 

= 1 — s jra" + ... — 1 — Ha n + ... , H positive. 

<p (a) n + 3 JJ r 

Put these in the differential equation in e of Art. 77 ; it 
gives 

B (m s + 5 m) a” -4 — GAffa*^ + ... = 0. 

Neither m nor B can equal zero, because then the second 
term of e only merges into the first. Nor can m = — 5, a nega- 
tive quantity. Hence the first term will not vanish of itself. 
But we may make the first and second vanish together by 
putting n = m and B ( m ’ + 5m) = GAIL Hence B must be 
positive. And therefore near the centre s increases towards 
the surface. 

In thus increasing, suppose it attains a maximum , and then 
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decreases. At this point ^ = 0 ; and the equation of Art. 77, 
already used, gives 

!?= {* - ^)} ? 0 

This corresponds to a minimum. Hence g does not attain a 
maximum, and therefore it continually increases from the 
centre to the surface. In the above we have assumed that 
<f> (a) is greater than pa 3 . This appears 

V <f> (a) = 3 1 ‘pW'da — pa 3 — f V d £ da, 

and ^ is negative by hypothesis. 

Prop. To calculate the numerical value of the ellipticity of 
the surface in the case of the Earth. 

80. In order to do this it is necessary to find the values of 
q& and tan q& at the surface. Let n be the ratio of the density 
of the surface to the mean density of the Earth. Now the 
mean density 

f 8 f s m 3(1 f a 1.) 

= J bar pa 1 da -f- j iva^da' = j— - cos ga + sin jaj- ; 

. 1 _- 3 (i ? a V 
n q V \ tan q&J ’ 

3nq& 

.’. tan q& = . 

1 3 n— q& 

If we take the mean density double of the density at the 
surface (see Art. 58), then 

tan qa - 3ga 

tan? a — 3 _ 2?v > 

which is satisfied by qa = 2'4576 = 140* 45'. Then 

tan ja = — 0*812, 0=4*0266, ffV = 6*0398, jV-s - 0 = 1*5. 
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Hence by the formula of last Article 

5m 1 m 

£ = Y 2*5266 = 1-01064 


1 

292 ’ 


since m = 


289 


If we take the mean density = 2-4225 x the density of the 
surface =2-4225 x 2 - 75 = 6'66, which accords very nearly with 
Mr Airy’s determination from the Harton Experiment (Phil. 
Trans, for 1856, p. 355, where it is 6-565), the equation for 
finding qa is 


tan qa = 


1 - 


qa. 

0-8075jV ’ 


which is satisfied by qa = 2‘618 = 150°. Then 


tan qa — — 0*57735, z = 5-5345, jV = 6-8539, 
2*a* -i- z = 1-2384, z + <fa’ = 0-8075. 


In this case 

5m 1‘4225 5m 1 m 
e = ~2 3-7729 = ~2~ 2*652 = 1-0608 


1 

306*6 ’ 


putting m = — 


289' 


81. In the course of the last eight Articles we have de- 
veloped the following conditions, which must be satisfied if 
the Earth has derived its present general form from being in 
a fluid state. (1) The direction of gravity must everywhere 
be perpendicular to the surface. (2) The form of the surface 
must be an oblate spheroid, with its axis coincident with the 
axis of revolution. (3) An additional test, though not abso- 
lutely infallible yet invested with a large degree of probability, 
is that furnished by the result of Art. 80, by assuming a law 
of density of the strata which is of itself h jyriori very pro- 
bably true, that the value of the ellipticity is not very dif- 
ferent from . We shall see in a future Chapter that the 
actual measurement of the form of the Earth by means of 
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trigonometrical operations meets all these conditions. It is 
found by measurements in widely separated countries, that an 
ellipse of the kind described can be drawn in the plane of the 
meridian of any place, cutting the plumb-lines at all the 
stations where it is examined at right angles ; and the ellip- 
ticity of this ellipse is almost exactly equal to There 

are local deviations from this law, arising from local causes, 
which are produced by the variations of the surface of the 
Earth and probably of the interior of the solid crust also. 
But the average line is this ellipse. Since the variations of 
the Earth’s surface, in mountains and valleys and extensive 
oceans, are palpable, and must have arisen since the Earth 
ceased to be fluid and assumed its general form, the fact that 
deviations from this ellipse are found in the level-curve while 
the average curve is still this ellipse, is rather confirmatory of 
the theory of original fluidity than otherwise. 

The probability of the truth of the law of density made 
use of in the previous calculations is strengthened by the 
value of Precession which it leads to. 


Prop. To test the laio of density used above by the amount 
of Precession of the Equinoxes wh ich it leads to. 


82. The Annual Precession 


C-A 

G 


3ra' r (, n " 2 

— cos id 1 + f 
n \ n 


1 - - sin* i 

“ 


180°, 


I = obliquity of the ecliptic = 23° 28' 18", i = inclination of 
Moon’s orbit to ecliptic = 5° 8' 50", n and n are the mean 
motions of the Earth round its axis and round the Sun, and 
their ratio = 365*26, n" the mean motion of the Moon round 
the Earth = 27*32 days, v = ratio of masses of Earth and 
Moon = 75. (See Mechanical Philosophy, Second Edition, 
Art. 470: also, changing the notation, Airy’s Tracts, Fourth 
Edition, p. 213, Arts. 36, 38.) Substituting the above quan- 
tities, 

C — A 

Annual Precession = 16225"*6 — 7 , — , 

c 
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where A and C are the principal moments of inertia of the 
mass, the latter about the axis of revolution. To find these 
let xyz be the co-ordinates to any element of the mass, rdco 
he the polar co-ordinates to the same. Then the mass of this 
element = — pr 2 dp. da dr, p, = cos 0. Also 


f + z 2 = r* {1 — (1 — f) cos* o>} 


= r 2 

"2 | 
[_3 + | 

^-(l-^)cos 2 <»] 

0 

£C* + 2 s = r 2 

"2 1 
_3 + 1 

[I - (! - /**) sin 2 ®! 

0 




The terms are here arranged as Laplace’s Functions. (See 
Art. 39, Ex. 4.) 

,'.C — A=[ f f p{(x* + f) — (?/*+z 3 )}dpdadr 
J -i J o Jo 


=J j J pr* [g - /**) - jg - (1 - f) cos’ coj dp. da dr. 
Now r = radius of any stratum = a jl + e g — /i’jj (Art. 75) ; 

!. f "' dr= lL<‘^r da 

= 11 / 35 { , + 5 ‘(s “'*’)}] * 

= // f ‘ + ^( I -/*’)}** 1 

= CT (a) + f (a) g - p.'j suppose ; 


P. A. 
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C-A=y/r (a) J J Q - s') {-/S+ (l- /i*) coa* to] dfidat, 
by Art. 26, 

= 7 rf (a) J Q-/**) ( J ~ V) ^ = ^ ( a )* 

o_r 

Also C= <r (a), neglecting the small term \[r (a). 

Now -fr(a.)=jp = 3 aS ^ ( a ) ( € — Tj ) 


= sin ? a (e - j) *, by Arts. 77, 78. 

And putting p = — sin qa, and integrating by parts, 

Ob 

tr (a) = [ pa* da = Q f a 3 sin qada 
J o ^ 0 

/ a 3 3a s . ,6a 6 . \ 

= $ — cos ^ sm q& + cos ya — — , sin ja j 

<?a* . L 6 / 6 \ qa \ 

= Y 8, “ « a f - ?i‘ - 1 1 - jvJ tST^f • 


Hence substituting z, as before, 



Substituting for q&, z, e and m their values, this is found 
to = 0-00313593. 

.*. Annual Precession = 16225"'6 x C’00313593 


= 50"*8. 


Tlie value generally assigned to the Precession, from obser- 
vation, i 3 50"- 1. The almost complete coincidence of the 
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result of the theory with thi3 observed value is a remarkable 
evidence in favour of the law of density we have adopted. 

83. Mr Hopkins has endeavoured to ascertain how far 
the interior of the Earth may at present be fluid, by calculating 
the value of the Precession upon the supposition of the mass 
being a spheroidal shell of heterogeneous matter, enclosing a 
heterogeneous fluid mass, consisting of strata increasing ac- 
cording to the law we have used. In three memoirs in the 
Philosophical Transactions of 1839, 1840, and 1842, he enters 
upon a complete investigation of this subject. We will give 
the evidence upon which he rests his conclusion that the crust 
is very thick. 

Prop. To trace the argument drawn from Precession to 
show that the crust is of considerable thickness. 

84. Mr Hopkins has deduced the following formula (in 
which we have changed the notation to suit the present 
treatise), 



where P is the precession of the equinoxes of a homogeneous 
spheroid of ellipticity e, which by calculation = 57" nearly 
if e = 5 ^; P is the precession of the heterogeneous shell, the 
outer and inner ellipticities being e and g: this = 50' ,- l by 
observation. 

The success of the calculation depends upon a remarkable 
result at which he has arrived, that the precession caused by 
the disturbing forces in a homogeneous shell filled with homo- 
geneous fluid, in which the ellipticities of the inner and outer 
surfaces are the same, is the same whatever the thickness of 
the shell. It is therefore the same for a spheroid solid to the 
centre. The formula above given is the relation of the 
amounts of precession in two shells, one heterogeneous and 
the other homogeneous ; and, as the thickness is the quantity 
sought, neither of these amounts could be calculated, and 
therefore the relation expressed in the above formula would 

6—2 
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be of no avail. But in consequence of the property that the 
precession of the shell, when it and the fluid are homogeneous, 
is the same as that of the spheroid, this difficulty is overcome ; 
and P can be calculated without knowing the thickness, and 
therefore P will be known. 


We have shown (Art. 79) that the strata decrease in ellip- 
ticity in passing downwards : hence g' — e is never negative, 
and the fraction on the right hand in the above formula is 
never negative, and is never so large as unity : let it = /S. 
Hence 


£ 

e 



or £ is less than - e ; 

8 


and therefore, because the ellipticity decreases in descending, 
the thickness must be greater than would correspond with an 
ellipticity of the inner surface of the shell equal to 7-8 ths of 
that of tne outer surface. 


If solidification took place solely from pressure, the surfaces 
of equal density would be surfaces of equal degrees of solidity. 
If we use the formula for finding e in Art. 78, and make 
ja = 150°, and the mean density = 2*4225 times the superficial 


density (the second of the values in Art. 80), then if e = -e in 

8 

3 

the formula of Art. 78, we have, after reduction, a = -a, or 


the thickness equal to one fourth of the radius, or 1000 miles. 
If a smaller ratio of densities is used than 2*4225, the thickness 
is greater. (Mr Hopkins shows also that a ratio a little larger 
than 3 makes the thickness 1-5 th of the radius: but this ratio 
is too large. The ratio generally used is about 2*2). 


But solidification depends upon temperature, as well as 
upon pressure. In his third memoir (Phil. Trans. 1842), 
Mr Hopkins shows that the isothermal surfaces increase in 
ellipticity in passing downwards. If temperature alone re- 
gulated the solidification, these surfaces would be the surfaces 
of equal solidity. But since both pressure and temperature 
have their effects, the ellipticities of tho surfaces of equal 
solidity must lie between those of the isothermal and the 
equi-dense surfaces. Hence the surface of equal solidity at 
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any depth will be more elliptic than the surface of equal 
density at that depth : and therefore the inner surface of the 

.7 

solid shell, of which the ellipticity is - e, must be at a depth 

8 i 

corresponding to a stratum of equal density of smaller ellip- 

7 

ticity than - e, that is, at a greater depth than 1000 miles. 

8 


In the above reasoning 0 has been neglected. If its value 
be used, it strengthens the argument for a greater thickness 
than 1000 miles. 


We may, therefore, safely conclude that 1000 miles is the 
least thickness of the solid crust. In the calculation it has 
been assumed that the transition from the solid shell to the 
fluid nucleus is abrupt. This will hardly be the case. The 
above result will therefore apply to the effective surface, lying 
near the really solid shell. But in consequence of the ten- 
dency, as shown above, of every cause being to prove that the 
crust is really thicker than 1000 miles, we may safely take 
this to be its least limit. 


85. Professors Hennessy and Haughton have both written 
upon this subject: see Phil. Trans. 1851, and Transactions of 
the Royal Irish Academy , 1852. The first makes the thickness 
be between 18 and 600 miles. But in his calculation he 
assumes that the shell is so rigid as to resist, without change 
of form, the internal pressure which arises from the inner 
surface ceasing to be one of fluid equilibrium : an assumption 
which cannot be considered admissible. Moreover he sup- 
poses that in cooling the outer shell will contract less than 
the fluid nucleus ; which can hardly be true. 

Mr Haughton’s investigation is simply a problem of den- 
sities, and determines nothing whatever regarding the ratio 
of the solid to the fluid parts of the Earth. (See Philosophical 
Magazine , Sept. 1860.) 

Prop. To show what influence the present aspect of the sur- 
face of the Earth has upon the argument for the thickness of 
the crust. 

* 
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86. The following considerations are sufficient to show 
that the crust of the Earth must at the present date be very 
thick. 



The above diagram represents a vertical meridian section 
of Hindostan, the Ocean, and the crust of the Earth, through 
0 , or Cape Comorin. AbcD is the average form of the 
mountain mass : AB = 140, BG = 230, Bb— Cc = 2*5 miles : 
mn = t, Ar = t' the thickness of the crust below m (any point 
on the table-land) and A : arc AM = a, area of AbmM = K ; 
G its centre of gravity, Gg vertical ; rg = k, perpendicular to 
Gg ; Mm — h ; h and k the middle points of mn and Ary he 
perpendicular to mn ; re —y, to he. 

The mass Mr is held in equilibrium by its weight, the 
downward pressure of the overlying mass MA, the upward 
pressure of the fluid below, and the force of adhesion at 
the joints mn and Ar. Since the crust has, by hypothesis, 
been formed by the solidification of the fluid, its density at 
any point will be very much the same as the fluid was at 
that point. We will at present assume it to be the same. 
Hence the weight of Mr = the upward pressure of the fluid, 
and the weight of Mb A tends to break the crust, and is sus- 
tained by the adhesion at the joints. Let G be the length 
of rock, of a unit section, the weight of which equals the 
average force of adhesion on a unit of surface. G — l-5th 
mile, may be considered to be the greatest limit of G (see 
Phil. Trans. 1855, p. 102). 

If the point m sink, the joints Ar and mn will open at A 
and n, ancl an opening will take place at some other point on 
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the left. The equation of moments of the forces acting on 
Mr taken about r is 


K.lc — C.t.y + C.t'. \t' 

= C .t{r — \t + h— {r — <') cos a} + %0. <'*; 

.*. 2K. & = {<'*+ 2tt' cos a — <* + 2 (r vers a + ^) <} (/...(l). 

Take the case of m being at c, then k — rg — 255 miles: 
a = 5° 19', cosa = 0’9957, rversa = 17‘2 miles: take K=Cb 
{omitting ABb), this = 230 x 2 - 5 = 575 square miles ; h — 2‘5 ; 

.*. <'* + 2 t!t - <* + 39'4< = 1466250 = (1212)* nearly. 

If t be very small, <'=1212 nearly; this is a condition 
which no law of cooling could, bring about. Also t' cannot 
be small, otherwise <* would be negative. If < = <', then each 
is greater than 800 miles. 

Formula (1) may be applied to find the least thickness of 
the crust beneath 8, any point on the Ocean south of Cape 
Comorin, to prevent its bed Os being broken up by the lava 
from below. Make 0 the centre of moments : Op = <', sq = t, 
OS = a, Ss = — h. Suppose that the depth of the Ocean in- 
creases uniformly with the distance from 0, and is 3 miles at 
25° distance from 0, i. e. between Madagascar and Australia ; 
then h = 14 sin ^a. Also K = area OSs ; and K. k, the moment 
about 0 of the several elementary portions, = 18’7r s sm s -|a by 
integration. Hence formula (1) becomes 

<'* + 2 «' cos a — f + 2 (r vers a — 14 sin Ja) < = 187^ sin® ^a. 

If a is so taken, that the coefficient of < may be neglected, 

<'* + 2<’< cos a — <“ = 187r*sin , ^a. 

As before, neither < nor <’ can be small. If < = <' ; then 
each = r sin V 93‘5 sin = 1000 miles, when a = only 10°. 
In this case the coefficient of < = 118, which may be neglected. 
Therefore, as before, the thickness must be very great to pre- 
vent the crust being broken through. 
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It has been assumed, that the density of the crust is every- 
where the same as the fluid from which it was formed by 
solidification. Suppose, however, that it is more dense, then 
the tendency of the crust to break, in the first case, will be 
greater than we have made it, though the tendency in the 
second will be less. The reverse will be the case if the crust 
is lighter than the fluid from which it is formed. So that 
in any case a fracture must take place somewhere, either 
beneath the mountains or beneath the ocean, unless the thick- 
ness is very great. As both the mountains and the ocean- 
bed retain their positions of elevation and depression, we have 
no alternative to choose but that the thickness of the crust is 
very great. 

87. The result of the whole proves that the crust must be 
very thick : and, as Mr Hopkins’s calculation appears to be 
free from objection, and in fact to be the only one on which 
any reliance can be placed, we may conclude that the thickness 
is at least 1000 miles. 

The present form of the surface in mountains, table-lands, 
continents, and oceans has been, no doubt, acquired from a 
process of expansion and contraction which the crust has 
undergone during the ages since it was first consolidated. 
Geology teaches us that these elevations and depressions of 
vast regions are at this present day going on. We may, 
therefore, fairly conclude — especially with this evidence that 
the crust is so thick — that the present varieties of the Earth’s 
contour have grown from this cause, and have not arisen in 
any way from the operation of hydrostatic principles. 

This does not in any way contravene the hypothesis that 
the Earth was once a fluid mass, and has received its general 
figure from that condition. The fact that its mean form, as 
measured by geodetic operations, coincides with the fluid-form 
calculated upon an assumed, but (h priori) very probable, law 
of density, is an unanswerable argument in favour of the hypo- 
thesis of original fluidity. And the coincidence of the calcu- 
lated amount of Precession, on this law of density, with its 
observed amount, is a very strong evidence that that law of 
density is the law of nature. 
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CHAPTER II. 

THE FIGURE OF THE EARTH, ON THE SOLE HYPOTHESIS 
OF THE SURFACE BEING A SURFACE OF EQUI- 
LIBRIUM AND NEARLY SPHERICAL. 

88. Upon the hypothesis of the Earth being a fluid mass 
it was shown by Clairaut, in his celebrated work Figure de la 
Terre, published in 1743, that the increase of gravity in pass- 
ing from the equator to the poles varies as the square ot the 
sine of the latitude, and that a certain relation must neces- 
sarily subsist between the ellipticity and the amount of 



simpler hypothesis of the surface only being a surface of equi- 
librium, and the interior being solid or fluid, but consisting 
of strata nearly spherical. Professor Stokes, in an investiga- 
tion published in the Cambridge Philosophical Transactions 
for 1849, has done the same, without making any assumption 
whatever regarding the constitution of the interior of the mass, 
but assuming only that the surface is a spheroid of equi- 
librium, of small ellipticity. The present Chapter is borrowed 
wholly from Mr Stokes’s investigation. Clairaut’s Theorem 
is valuable as it gives us the means of determining the ellip- 
ticity by means of pendulum oscillations, the times of which 
measure the force of gravity at the several stations where 
experiments are made. 

Prop. To find the law of gravity at the surface of a sphe- 
roid of equilibrium, and of small ellipticity. 


* 
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89. Let V be the potential of the mass. Then because 
the surface is a surface of equilibrium, 

const. = V+ \v? (1 — (j?) r*. 

% 

By a process precisely like that in Art. 23 we have, for an 
external point, 

d*.rV , d t . dV\ 1 d*V n 

r df ^ dp) + l-iSda'~ 0 ' 

Let Fbe expanded in a series of Laplace’s Functions, 
i^)+ v t + ••• + Fjd* ••• „ 


Then since the above equation is linear with respect to V, 
and a series of Laplace’s Functions cannot equal zero unless 
the Functions are separately zero (see Art. 35), we have, by 
substituting the above series for V and remembering the con- 
dition given by Laplace’s Equation, 

fP r v, 

rf 4P-*(* ,+1 ) F ‘ =0; 


w 

v^fi+zy, 


where W { and Z t are independent of r. The complete value 
of V becomes 


WWW 

V= + -^r + -^r + ... +Z B + rZ t + r i Z a + ... 

Now V evidently vanishes, from its very definition, when r 
is infinite. Hence Z a = 0, Z i — 0, Z a = 0 . . . 

WWW 
... F=~° + -^r‘ + ^*+... 
r r r 


If the surface were spherical, this expression would be re- 
duced to its first term. Hence in our case W t , W t ... must 
be all small quantities of the first and higher orders. 
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Substituting this in the equation of equilibrium and putting 
r = a{l + « (!—/)}, 

we have 

const. = ^°|l - e Q — M a )| + 5 r + -Jr + •••+ \ ^ (1 -/**)• 

Equate the sums of Laplace’s Functions of the same order 
to zero ; 

TF 0 = a const. — 5 w?*a 8 , W x = 0, 

O 


K = (*>-!««)»' (1 -s). 


w,=o, w,=o... 

■■■ r-Z+(Ws-iM)$(l-y). 

Let g be gravity. Then since the angle between the 
radius vector r and the normal varies as the ellipticity and 
therefore its cosine must be taken = 1, the value of gravity 
dV 

is — —j- — the part of the centrifugal force resolved along r 




3 ( ' - 1 ««) p (5 (5 + 5 " '**) • 


Substituting for r and omitting small quantities of the second 
order, 




W,e-iwW 

a 2 


(r"') • 


The first portion of this is evidently the mean value of g, 
because if g be multiplied by an element of the surface and 
integrated throughout, the latter part will disappear. Let 
the first part be O ; therefore also w* . a = m . G, since m is 
small; 
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+ («-§*)(!-,.•)• 

= O |l — m— ejj jl + sin* (latitude) j . 

Hence the increase of gravity in passing from the equator 
to the poles varies as the square of the sine of the latitude ; 
and also 


polar gravity — equatorial gravity 
equatorial gravity 


+ ellipticity 



= - x ratio of centrifugal force at equator to gravity. 


90. This is Clair nut's Theorem; which is thus demon- 
strated by Professor Stokes without making any assumption 
regarding the interior of the mass. Nothing can be inferred, 
therefore, from any numerical value we may obtain for gravity, 
and therefore for the ellipticity, by pendulum experiments, 
regarding the Earth’s mass having been originally fluid or 
not. 

For a valuable and interesting account of pendulum ex- 
periments made in places in all latitudes, and the result re- 
garding the Figure of the Earth, we must refer the reader to 
Major General Sabine’s work on the subject, Account of Ex- 
periments to determine the Figure of the Earth by means of the 
Pendulum vibrating Seconds in different Latitudes. London, 
1825*. The ellipticity thus deduced is sis', rather greater 
than that obtained, by the geodetic and other methods. In 
consequence of the irregularities of the surface of the Earth 
the experiments with pendulums need various corrections 
before they can be properly applied to determine the ellip- 
ticity. The principal ones depend upon the elevation of the 
station above the sea-level (for which Dr Young gave a formula 
of correction, see Phil. Trans, for 1819), and the excess or defect 

* See also his latest remarks on the subject in the Notes to his translation 
of the Comm, Vol. JTV. Part I. 
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of matter in table-lands or tbe sea in the neighbourhood of the 
station. If the station, for instance, be on a rock in an island in 
the midst of a sea, such as St Helena, the correction for this 
second disturbing cause will be different from what it would 
be for a station at the same elevation from the sea-level in 
the midst of a continent. This effect depends, as may be 
gathered from Art. 52, not upon the height of the station 
above the sea-level, but upon the excess or defect of matter 
however arranged. Professor Stokes has fully considered the 
influence of these causes of derangement in his Paper above 
referred to. He shows, that the effect of these corrections for 
the irregularities of the surface, and for the different elevations 
and other local circumstances of the stations where the experi- 
ments are made, is to reduce the value of the ellipticity, and 
make it nearer to 

91. Cor. 1 . The investigation in the last Proposition 
gives, after making substitutions, 



E being the mass of the Earth. W a is evidently equal to the 
mass, because as r becomes infinitely great the second term 
vanishes with reference to the first, and we know that in that 
case the value of the potential must be the mass divided by 
the distance. 

92. Cor. 2. Laplace first pointed out that the ellipticity 
of the Earth would have an effect upon the Moon’s motion. 
The expression in the last Article leads to the following for- 
mula for the change in the Moon’s latitude produced by this 
cause, 

■2^( 6_ l)“ in2/3in(,1< + e) ' 

where a = distance of Moon, h = mean motion of the node, 
n — mean motion, e = the epoch, or longitude when t = 0, 
I the obliquity of the ecliptic. (See Mechanical Philosophy , 
Second Edition, Art. 556 ; also, changing the notation, Airy's 
Tracts , Fourth Edition, p. 188, Art. 84.) This has been 
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shown by Burg and others to equal — 8" sin (nt + e ) ; also 
h = 0"0040217«. Hence after all reductions 


€ _ 1*0-0015474, f = 5^8 = 0-0017476; 


e = 0-0032950 = r nearly. 

303 J 


This method, as was the case with the pendulum experi- 
ments, only determines the ellipticity, but gives no evidence 
on the subject of the original fluidity of the Earth. The near 
agreement is remarkable. 

93. The spheroidal form of the Earth’s surface and the 
circumstance of its being a surface of equilibrium afford us 
more information, as Professor Stokes has shown, regarding 
the distribution of matter in the interior of the Earth’s mass. 


Prop. To show that the centre of the Earth's mass coincides 
with the centre of its volume, and that the axis of rotation is one 
of the principal axes of the mass, as a consequence of the form 
of the surface being a spheroid of equilibrium. 

94. By Art. 18 we have the potential of the mass with 
reference to an external point {fgh), or 

pdxdydz 

K/- *)* +{ff- yT + ( A - *)*]* ’ 

the integrals extending throughout the whole interior of the 
Earth. Let pdxdydz=dm, and X, p, v the direction-cosines of r 
the radius vector to (fgh) ; then f— \r, g = pr, h = vr\ and 
expanding the radical according to the inverse powers of r, 
we get 



V= \ fjj drn + p 2.X j' jjxdm + ^2 (3\ 2 - 1 ) [ jjx'dm 
+ p % . Xp IJj xydm + ..., 


S denoting the sum of the three expressions necessary to form 
a symmetrical function. 
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Comparing this with the value of Vin Art. 91, 

I ffdm = E, and also 

jjjxdm = 0, j jjydrn = 9, J fj zdm = 0 (1), 

i S . (3A 2 — 1) jj Jx * dm + 3S . Xfifj fxydm 

= ( e ~l m ) ( 2 )> 

and other equations from the succeeding terms. 


Equations (1) show that the centre of gravity of the mass 
is at the centre of the spheroid of revolution, or the centre of 
the volume. 

With regard to equation (2) we may observe that Xfiv are 
tied by the relation X 2 + /u? + v 1 = 1. If then we insert this for 
1 in the equation, so as to make it symmetrical with regard to 
X, fi, and v, and equate the several coefficients, we shall 
obtain 


S . (x 2 — i // — i j/ 1 j |"|j oddm + 3% .Xfi J ffxy dm 
which leads to 


f f fxydm — 0, fffyzdm = 0, fffxzdm = 0 (3), 

fjjx*dm-ljfjy a dm-lffjz*dm= 

lj!y‘dm-lfjf* t dm-lfjfa?dm= f Ea? 

jfjz'dm -Iff fx'dm - \ fffy'dm =- § (e ■ - 5*#) Ez* 


L ...(4). 
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Equations (3) show that the co-ordinate axes are principal 
axes ; and that therefore the axis of revolution is a principal 
axis. 

The last three equations (4) show, first, that 

JJJ x*dm — JJJy*dm, 

or, as might be anticipated, the moments of inertia about the 
axes of x and y, in the plane of the equator, are the same. 
They give also if, as before, A and C are the principal 
moments, 

C — A = JJJ r j?dm —JJJ z *dm = ^ ~ -Eh 1 . 

Were we able to deduce C also, without making any hypo- 
thesis regarding the internal condition of the earth, the Pre- 
cession might be obtained under the same circumstances. In 
that case, the agreement between the amount of precession 
calculated upon the assumed law of internal density and the 
observed value (pointed out in Art. 82) would furnish no 
proof of the truth of that law. But as the precession cannot 
be thus independently calculated, the result in Art. 82 affords 
a strong argument for the correctness of the law of internal 
density which we have adopted. 
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CHAPTER III. 

THE FIGURE OF THE EARTH, DETERMINED BY GEODETIC 
OPERATIONS. 

95. In a previous Chapter we have shown that if the 
Earth he considered a fluid mass the form of the surface will 
be an oblate spheroid of small ellipticity, its axis coinciding 
with the axis of revolution, and the surface being everywhere 
at right angles to the direction of gravity ; and further, that 
upon assuming that the density of the strata varies according 
to a certain very probable law, the ellipticity = nearly. 

In this Chapter we propose to submit this to the test of 
measurement, by inquiring whether an ellipse can be found 
with its axis coinciding with the axis of the Earth and cutting 
the plumb-line at stations along it at right angles ; and whether 
the ellipticity of that ellipse is ~ . 

The method of doing this is as follows. A base-line, about 
5 or 6 miles in length, is measured with extreme accuracy, near 
the meridian, the curvature of which we are to find. By a 
series of triangles thi3 base is connected with a number of 
stations in succession lying near the meridian, the angles and 
sides of which are calculated or observed, as the case may be. 
Thus a connexion is established between the original base and 
a second base at the termination of the chain of triangles, 
and the length of this second base obtained by calculation. 
It is then measured, as the first was, and by a comparison of 
the calculated and measured results the correctness or not of 
the operations is tested. This having been satisfactorily per- 
formed, the projections of the sides of the triangles on the 

P. A. *7 
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meridian are found, and their sum gives the length of the 
meridian arc between its two extremities. The latitudes of 
these extremities are then observed with great care, and from 
these data the form of the ellipse, of which the arc is a part, is 
found by the principles of conic sections, as we shall now 
show. 


Prop. To find the length of an arc of meridian in terms of 
the amplitude, the semi-axis major, the ellipticity ( the ellipticity 
being small), and the middle latitude. 

96. Let l and l be the latitudes of the extremities of the 
arc, m the mean of these or the middle latitude ; X the ampli- 
tude of the arc or the difference between the latitudes ; a, b, and 
e the semi-axes and ellipticity ; s the length of the arc, r the 
radius vector, and # the angle r makes with the major axis. 
Then 


1 cos*# , sin*# , , a* 

^5 — I -gr~ j t an l — i ; 


|jtan #; 


1 a* cos* l + V sin* l . . 

•P“7S?7+JWJ> P“tt“g J = a(l-<), 

r = a (1 — e sin*Z), neglecting e’... 

— — — 2ae sin l cos l, ^ = 1 — 2e + 4e sin* l : 
dl dl 


ds / , rf#“ dr ‘ . . s A 

jrv^ + ^ =a(1_2e+3esmZ) 

= o( 1 -^e-|ec°s2z); 

/. s = a |^1 - iej (l — T) — |e(sin2Z— sin2Z')j 
1 3 

= - (a + l)\ — - (a — b) sinXcos 2m. 
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97. Cor. If X be small, not exceeding 12°, we may put 
sin X = X in this formula ; then 


8 

X 


a + b 
2 


a — b 


cos 2 m. 


Prop. To obtain formulae for finding the semi-axes and 
ellipticitg, when the lengths, amplitudes, and middle latitudes 
of two small arcs are known; and to ascertain what arcs are 
adapted to give the best results. 

98. Let sXm, s'X'm be the lengths, amplitudes, and mid- 
dle latitudes ; 


s a + b 
*** X- 2 



cos 2 m, 


V 


o + 5 
2 " 


„ o — 5 . , 

3 — — cos 2m ; 


a — b 
•** ~ 2 ~ 


1 — — 

1 X~X' 

3 cos 2m — cos 2m ’ 


s s 

, , - cos 2m — --;Cos 2m 

a + b V 

2 cos 2m' — cos 2m ’ 


by which a and b and therefore e are found. 

The effect on the axes of any error in the amplitudes will 
be found by differentiating the above formulae. In the deno- 
minators of the resulting expressions the quantity 

cos 2m — cos 2m' 


will appear. The errors in the axes consequent on errors in 
the observed amplitudes will, therefore, be least when this 
quantity is a maximum. Suppose one arc is chosen ih the 
southern half of the quadrant, cos 2m is positive ; then 

2m'=180°orm' = 90° 

will give the best result. Suppose one arc is in the northern 
half, cos 2m is negative ; then 2m' = 0 will give the best 
result. Hence the nearer one arc is to the pole and the other 
to the equator, the less will errors in the data affect the calcu- 
lated form of the ellipse. This will be illustrated in the fol- 
lowing examples. 

7—2 
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99. Ex. 1. Compare the two parts of the Indian Arc 
from Kaliana (lat. 29°. 30'. 48") to Kalian pur (24°.. 7 . 11"), the 
length being 1961157 feet, and that between Kalianpur and 
Damargida (18°. 3'. 15"), the length being 2202926 feet. 

X = 5*. 23'. 37"= 19417", V = 6°. 3'. 56" = 21836", 

2 m = 53°. 37'. 59", 2 m = 42*. 10'. 26", 

.*. -j (a - 5) mx 54456, \ (a + 5) = 20929789 feet, 

** 2 

a = 20984245, 5 = 20875333, e = — -. 

1 JO 

Ex. 2. Compare the two parts of the English Arc ; viz. 
from Saxaford (60°. 49'. 39") to Clifton (53°. 27'. 30"), measuring 
2692754 feet, and from Clifton to Southampton (50°. 54'. 47"), 
measuring 928774 feet. 

X = 7°. 22'. 9" = 26529", X' = 2°. 32'. 43" = 9163", 

2m = 114°. 17'. 9", 2m = 104*. 22'. 17", 

/. i (o - b) = 59419, J (a + 5).= 20863630, 

J 2 

a = 20923049, 5 = 20804211, e = J~. 

’ 176 

Ex. 3. Compare the arc between Kalianpur and Damar- 
gida with that between Clifton and Southampton. 

X = 6*. 3'. 56", X' = 2®. 32'. 43", 

2m = 42*. 10'. 26", 2m' = 104°. 22'. 17", 

.*. I (a - 5) = 33094, \ (a + 5) = 20882770, 

a= 20915864, 5=20849676, e = ^. 

olo 

It will be seen in these examples that when the arcs com- 
pared are near each other the resulting ellipticity differs much 
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from that deduced by the fluid theory : but when they are 
more distant from each other, as in the third example, the 
result is far more accordant. If there were no errors in the 
data, viz. in the observed amplitudes and measured arcs, the 
results ought to come out in complete accordance with each 
other, if the figure of the Earth be truly spheroidal ; for the 
formulae are sufficiently exact for this purpose. 

Prop. To explain the cause of the ellipses, determined from 
the several pairs of arcs, differing from each other. 

100. We have assumed, (1) that the meridian arc is an 
ellipse, that being the form which it would have were the 
Earth fluid : (2) that the plumb-line at all stations of the 
meridian is a normal to this ellipse. These suggest in what 
direction we are to look for an explanation of the discrepancies 
in the results. 

First. It is obvious that the form of equilibrium no longer 
actually exists, as all the variety of hill and dale, mountain 
and table-land and ocean-surface, sufficiently testifies. Geology 
teaches the same more generally and philosophically. Exten- 
sive portions now dry land were once at the bottom of the 
ocean, receiving the fossil deposits and burying them in the 
detritus of rocks, which time wore down, to become, as they 
are now, the records of their own history. Changes of level 
must therefore have taken place on a large scale. Land- 
marks in Scandinavia, the temple of Serapis at Puzzuoli, the 
ancient and recent coral-reefs in the Pacific, as pointed out 
by Mr Darwin, all testify that these changes of level are still 
slowly going on. It has been suggested, with great pro- 
bability, to be caused by the expansion and contraction of 
vast portions of rock in the interior of the Earth arising from 
variations in temperature produced hy chemical changes. 
Whatever the cause, the fact is certain. The Earth’s form 
can no longer be a form of fluid-equilibrium, although the 
average form is so. 

Secondly. The plumb-line may not in all cases be per- 
pendicular even to the mean ellipse. Local attraction is suf- 
ficient to produce material errors in the vertical, and therefore 
in the amplitudes determined by meridian zenith distances 
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of stars. For instance (Art. 56), an error as great as 5" was 
discovered at Takal K’hera in Central India by Colonel 
Everest, arising from the attraction of a distant table-land. 
Sir Henry James has shown that a deflection of about the 
same amount occurs at Arthur’s Seat, Edinburgh (Phil. Trans. 
1857). We have mentioned that the attraction of the Him- 
malaya Mountains produces a deflection amounting to as much 
as 28" at the northern extremity of the Great Indian Arc 
(Art. 61). We have calculated elsewhere (see Art. 62 and 
Phil. Tram, for 1859) that the deficiency of matter in the vast 
ocean south of India causes such deflections as 6", 9", 10"‘5, 
19"‘7 at various stations: and (Art. 63), shown that it is not 
improbable that extensive but slight variations of density 
prevail in the interior of the Earth, .the causes of which are 
not visible to us as mountain masses and vast oceans are, 
sufficient to produce errors in the plumb-line quite as great 
as and even greater than most of those already enumerated. 
These seem abundantly to account for the variety in the cal- 
culated semi-axes and ellipticities in the last Article, derived 
as they are from uncorrected observations. 

101. Mr Airy has entered very thoroughly into a com- 
parison (see Figure of the Earth, Encyc. Metrop.) of the various 
arcs measured in different parts of the world. He has used 
them according to their importance and value, as determined 
by the circumstances under which they were measured and 
observed. His result satisfactorily shows that the ellipticity 
of the mean spheroid is about ^ . The conditions, therefore, 
required for supposing the Earth to have received its present 
average form from having been once in a fluid state, are 
altogether satisfactorily fulfilled. 

The same result has been obtained by another process, first 
used by the late M. Bessel and adopted by Captain A. Clarke, 
R. E. in the Volume of the British Ordnance Survey. This 
method we shall now explain, first introducing one or two 
propositions which we shall require for its application. 

L’et the form of the meridian line be such, that 
p = A + 2B cos 2l + 2 C cos 4Z 
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is the radius of curvature at a point of which the latitude 
is Z. 


Prop. To prove that if the meridian he an ellipse, 

6CA = 5B\ 

102. Let p'x'y he the radius of curvature and co-ordinates 
to a point in latitude l, in an ellipse, 


a? y* s b* 

a* + b i = 1 ’ 6 =1 ~a" £ = 


tsnZ= -$> />'-( 

From these we obtain 


1 + 


*4 


! . d *y 


dxV - dx'* ' 


x' — a cos l (1 — e s sin* Z) _i , 


y' = a sin l (1 — e*) (1 — e* sin* Z)~ l , 
p' = a (1 — e*) (1 — e* sin* Z)“*. 


Expanding this last, neglecting e*.... 

p = a (1 — e s ) (l+|e a sin s Z + ^eWZ) 

A O 

= “{ 1 -; e ’-| e ‘-(! e ’ + 4) C0s2!+ s e ‘ COB4 f 

Comparing this with p' = A + 2B cos 2Z 4 2 G cos 4 Z, 

••• 6 ^= 4 4 a ’4T JJ, = 5S '- 
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103. COR. The expansions of x and y are as follows : — 

*'= ° {(* ■ + V + l e ‘) “ s l ~ G e ’+ xi5 e ‘) cosW+ jf/ 008 , 

y=a|( 1 -|e 1 -^) 8 m!-(I«'-ie') 8 m3(+ie‘8in5!-. 

Prop. To find an expression for measuring the departure 
of the curve of the meridian from an ellipse at any point , when 
the meridian is not elliptical. 

104. Let xy be the co-ordinates to any point of which l is 
the latitude, and the radius of curvature as above. 

-vr dy dfy , 7 dl 

Now cot l = — -f - , .*. = cosec l ~r ; 

dx dor dx 

^ = — p sin l, and ^ = p cos l; 
x = (A — B) cos l + ^ (2? — C) cos 3 1 + \c cos 5l, 

y = (A + B) sin l + ^ (B + C) sin %l + \c sin 5 1. 

o o 

Let o and h be the semi-axes of the curve, whether it be 
an ellipse or not. Hence these values of x and y give 

a-AS + \(B-0+\c=A-\B-y, 


b = A + B-\{B + C)+\c=A + \B-?- b C, 

e* = 1 - ^ ^ ( l +| 2 )» neglecting B.C &c. 

If we put these for a and e % in the expressions for x and y 
in Art. 103, we shall have the co-ordinates of a point (lati- 
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tude l) in the ellipse constructed upon the same axes as those 
of the actual curve. After reduction, we get 


X — X 


y -y' - (°- ! j) (A sin 1 + s sin il + 1 8in M ) • 

Let hs and hr be the distances between the points (xy) and 
(• xy ) measured along the arc of the ellipse and the normal ; 

5s =— [x— x') sin l+{y—y') cosZ=-^ ^C— ^ ^-) sin4£, 
hr = (x — x ) cos ? + (y — y') sin l = ^ ^(7— ~ sin 8 2 1. 


Prop. To obtain a formula for correcting the amplitude of 
an arc , so as to make its measured length accord with a given 
curve. 

105. Let s be the length of the arc and p the radius of 
curvature as before ; then, by integration, 

s = Al + B sin 2l + ^ C sin 4 1 + constant. 

Let l — \d>, l+\(f> be the limits of s, l being the latitude 
of the middle point, and cj) the amplitude of the arc ; 

s = Af> + 2Bcos2l sin <f> + (7cos4Zsin 2 (f>. 

Suppose now that x x x{ are the small corrections which must 
be applied to the observed latitudes, l — \(f>, l + \<f>, to make 
them accord with the measured length s. Then l— \<f> + x l 
and l + %<f> + x' and <£ + #,' — must be put instead of 
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l — \$, l + and <f> in the above formula. Hence, neg- 
lecting the squares of small quantities, 

s = A (<f> + — » t ) + 2i? cos 2 1 {sin <f> + («,' — x t ) cos <f>] 

+ C cos 4 1 {sin 2<f> + 2 (a>/ — aj cos 2<f >] ; 

.*. (x’ — x t ) (A + 2 B cos </> cos 2l) 

= s — A(j> — 2B cos 2? sin <\>—G cos 4Z sin 2<j>. 


Put A + 2B cos <f) cos2l = A +fj.; 



sin <f> cos 2 1 — sin 2<f> cos 4?. 


Let 


1 1 

A ~ 2089000 



-iLl 

10000 / 


2 B = IF C = Z 
A ~ 200 + 10000 ’ A ~ 10000 ’ 

and the above formula will become 


® l , = 7w + aU+(3V+ 7 Z+ x v 

where m, a, /3, 7 are functions of the observed latitudes, the 
measured length, and numerical quantities only. 


106. As an example which the student may work out for 
himself, the following is selected from the Volume of the 
Ordnance Survey. 


Station. 

Observed 

Latitudes. | Amplitudes. 

Measured 
Arcs in feet. 

Damargida 

18°3' 15"*292 



Kalianpur 

24 7 11 *262 

6° 3'55"*970 

2202904*7 

Kaliana 

29 30 48 *322 

11 27 33 *030 

4164042*7 
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If x t be the correction for Damargida, then the formula, 
when the numbers are substituted, will give these correc- 
tions : 

For Kalianpur...— 4-063 + 2'1831 U+ 1-6212 V+ 0-4285Z+ x , , 

... Kaliana + 0-365 + 4*1251 U+ 2*7741 F- 0-7213 Z+ x v 

Also Br = - 29-01 - 1'16 F- 557-07Zfeet. 

Prop. To explain the process by which the mean figure of 
the earth is obtained from the observed latitudes and the mea- 
sured arcs by the Principle of Least Squares. 

107. Suppose that we have a number of equations with 
numerical coefficients connecting a number of unknown quan- 
tities, less in number than the equations ; if the equations are 
true, the same values should come out whichever of the equa- 
tions we use in the process of elimination. In Physical 
Science it often happens that we have a problem of this sort, 
in which the numerical coefficients, being obtained from ob- 
servations, are not exact, but only approximate. If the cor- 
rect values of the unknown quantities were substituted, they 
would not exactly satisfy the equations, but small residuary 
errors will appear, differing according to the set of equations 
we select for elimination. It would seem, therefore, difficult 
to determine which of all the results actually obtained is 
nearest the truth. The late Professor Gauss discovered the 
Principle of Least Squares, which is of eminent service in 
such cases of perplexity. The principle is this ; that those 
values of the unknown quantities are nearest the truth which 
make the sum of the squares of the errors the least possible. 

In using this principle the Differential Calculus will evidently 
furnish us with exactly as many equations as there are un- < 
known quantities ; and the problem will be solved, with the 
nearest approximation to the truth attainable. 

In this manner the Mean Figure of the Earth may be de- 
termined. In the Volume of the Ordnance Survey eight arcs 
in Europe and India, consisting of 66 subordinate portions, 
have been used. In each arc the errors in the latitudes of 
the principal stations, which divide it into its subordinate 
portions, are calculated, as in the last Article, in terms of 
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UVZ and the unknown error (as *,) in one of the terminal 
stations, or in any one of the stations chosen as a starting 
point. The eight arcs will thus furnish 66 formulae of cor- 
rection, (similar to that for x[ in Article 105), involving 
eleven unknown quantities U, V, Z, x 1 ,x t ... x s . The nearest 
values of these are obtained by differentiating the sum of the 
squares of these corrections with respect to these unknowns 
and equating the results to zero. The process of calculation 
is very laborious. 


108. In the Ordnance Volume U= — 0*6937, F=1'4838, 
.£=0*3739; and these make x 

a = 20927197 feet, b = 20855493, e=-^-, 

M J U 


and the value of Sr in Art. 104 becomes 117*5 sin’ 2 1 feet, 
which shows that the greatest departure from the elliptic 
form is in latitude 45°, and equals 117 *5 feet. The correction 
of the latitude of Damargida (i.e. the value of x x in Art. 106) 
is — 0"*246, and the consequent corrections for Kalianpur and 
Kaliana are — 3"*578 and 1"*643, for the above mean values of 
U, V and Z. The above measures determine that curve which 
is nearest to the meridian, of all the curves represented by 
the general formula in Art. 101. It appears to be very 
nearly elliptical, bulging out but slightly in the middle 
latitudes. 


Prop. To explain the process of finding the Ellipse most 
nearly representing the observations. 

109. The process is precisely similar to that explained 
in the last Articles, C being first made equal to 5Zr-r6A, 
that the curve may be an ellipse. By Art. 105, 

z _ 5 / J . v V . z= _ ±__Z 

10000 24 V200 10000/ ’ “ 96 480* 

There will be only ten unknown quantities in this case. 
The corrections for Aalianpur and Kaliana in terms of x 1 , 
that for Damargida, are 


\ 
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- 4"-085 + 2"'l 831 TJ+ l''-6203 V+ x t , 

0-403+ 4-1251 U+ 2’7756 V+x 1 . 

The calculation makes Z7= — 0-3856, F= 1-0620. The 
value of Sr is zero when Z is substituted, as of course it 

should be. Also a = 20926348, b = 20855233, e = . Sir 

294 

Henry J ames takes as the final result for the mean figure of 
the Earth (see his Preface), 

a = 20926500, 5 = 20855400, e = — -. 

The corrections for latitude (in the example we have taken 
all along) are * 1 = 0"-050, for Kalianpur — 3''-156, and for 
Kaliana l"-810. These are the quantities by which, according 
to the Principle of Least Squares, the observed latitudes must 
be altered to make the measured arcs accord with the mean 
ellipse above determined. 


110. What has gone before leads to the determination of 
only the Mean Figure of the Earth. Any one meridian may 
possibly differ from this mean form owing to local causes, such 
as the rising or sinking of the surface from internal expansion 
and contraction of the materials of the crust, which may have 
taken place since the form ceased to be regulated by the laws of 
fluid equilibrium. Indeed the average result even seems to point 
out that some such change has occurred. For it appears in 
Art. 108, that the calculations in the Ordnance Survey Volume 
show, that there is a slight protuberance in the middle lati- 
tudes, even in the mean figure of the earth. 

If there be any local deviations from the mean figure fur- 
ther than this, India seems to present phenomena which 
would suggest, that these deviations must exist there if any- 
where ; and in that country an extensive and well-executed 
Survey has been carried on, which supplies us with data. 
The particular case of the Indian Arc has been used for illus- 
tration in the preceding references to the Ordnance Volume, 
because the formulae will now be of use in the following cal- 
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dilations. In India there are visible sources of error in the 

5 osition of the plumb-line which ought not to be overlooked. 

'he mountain-mass on the North, and the ocean on the South 
by its deficiency of matter, both tend to give the plumb-line a 
deviation northward and through different angles. The author 
has approximated to the amount of these deflections, as before 
stated in this Treatise. The following are the results he has 
obtained. (See Phil. Trans. 1859.) 

Deflections at Damargida, Kalianpur and Kaliana are : — 

Caused by the Mountains 6"*79 12"‘05 27"‘98 
Ocean 10 *44 9 *00 6 *18 

Totals 17 *23 21 -05 34 '16 

Errors in the amplitudes ... 3"‘82 ... 13"*11. 

If these be applied as corrections to the amplitudes in 
Ex. 1 of Art. 99, we have A = 5°. 23'. 37"+ 13''= 5°. 23'. 50", 
\'=6°.3'.56" + 4'' = 6°. 4'.0’’, and the formula! of Art. 98 
will give 

a =20906792, & = 20843795, e = — , 

004 

which is nearer the mean ellipse than the uncorrected data in 
Ex. 1 make it. 


111. Captain Clarke has suggested the following course. 
By the principle of least squares he finds the ellipse which 
differs least from the mean ellipse in form, and gives deflec- 
tions of tlje normal from the normal of the mean ellipse most 
in accordance with the calculated deflections. This he has 
done, taking account of mountain attraction only; the effect 
of the ocean on the plumb-line had not then been estimated. 
We propose now to go through his calculation, taking account 
of both these visible causes of disturbance. 

Prop. To determine the ellipse which most nearly accords 
with the mean ellipse in form, and at the same time most 
nearly meets the anomalies in India arising from mountain 
and ocean attraction. 
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112. Let Z,, Z s , Z s be the latitudes of Damargida, Kalian- 
pur, and Kaliana, corrected to the mean ellipse, so that (see 
Art. 109) the observed latitudes are Z, — 0"'05, Z a 4- 3"T6, and 
Z s — 1"*81. If, now, taking a general case, Z, + e,, Z a 4- e v Z a +e, 
are the latitudes for the three places referred to any other 
ellipse, thene, + 0 - 05, e a — 3*16, e s — 1*8 1 are the corrections 
which must be added to the observed latitudes to make them 
accord with the new ellipse. Hence by Art. 109, 

e, + 0*05 = se„ 

e a - 3-16 = - 4-085 + 2-1831 U+ 1-6203 V+ x x , 
e s + 1-81 = 0-403 + 4-1251 U+ 2‘7756 V+ x x . 

U= - 0-3856 + 1-85006, - 4-4446e a + 2-5946e s , 

V = 1-0620 - 3T098e, + 6-6056e a - 3-4958e a . 

Suppose that d v <Z a , d s are the angles of deflection caused 
by the mountains and the ocean. Then the ellipse which will 
most nearly satisfy the Indian Arc is that which makes 

(6,- < Z,)*+(6 a - < Z a )*+( ei - < Z I )* 

+ (1-85006, - 4-4446e a + 2-59466,)* 

+ (-3-1098e,+6-6056e a -3-4958e # ) s 

a minimum. By differentiation, with respect to e,, e a , e s , we 
obtain three equations, which after transformation become 

e, = 0-82493 d 1 + 0-30087 d % - 0-12583 cZ a , 

e a = 0-30086«Z, + 0-34199<Z a + 0-357 16^ 3 , 

6, = - 0-12584(Z, + 0-35715<Z a + 0-76873i a , 

and from these we find 

U= - 0-3856 - 0T3760J, - 0*03671<Z a + 0-17432 
V = 1-0620 - 0‘13808cZ, + 0’07484cZ a + 0-06325cZ a . 

The values of <Z,, <? a , d a are now to be substituted: they 
are 17"-23, 2l"-05, 34"*16 ; and they make 6, = 16"-25, 
e a =24"'58, 6 S = 31"*61 ; also U = 2*4255, F= 2-4189. 
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Hence the errors in the observed latitudes as affected by 
deflection, or e^O'^OS, e 3 — 3"*16, e 3 +l"*81, are 16"’30, 
2l"*42, 33''*42. These are less than the deflections by the 
small quantities 0"‘93, — 0"*37, 0"‘74. 


The values of U and V give the following results : — 


a = 20919988, 5 = 20846981, e 


1 

287 * 


113. We have thus obtained three different measures of 
the arc in question : viz. I. That derived from a comparison 
of the two portions of the arc together, the amplitudes not 
being corrected for local attraction ; II. The same comparison 
after the amplitudes are corrected for mountain and ocean 
attraction ; III. The ellipse obtained by least squares, which 
departs least from the mean ellipse in form and at the same 
time gives deviations of its normals from the normals of the 
mean ellipse as nearly as possible equal to the calculated 


deflections arising from local 
here gathered together : — 

a 

attraction. 

5 

The results are 
€ 

Mean Ellipse 

20926500 feet, 

20855400, 

(Art. 109), 

Arc I. 

20984245 „ 

20875333, 

4, -(Art. 99), 

Arc II. 

20906792 „ 

20843795, 

sb- (Art - uo) ’ 

Arc III. 

20919988 „ 

20846981, 

(Art. 112). 


Let 3a and 85 be the excess of a and b for each of these 
ellipses compared with the mean. 

Arc i. Arc n. Arc m. 

Hence 8a = 10’93 miles, —3*73, — T23, 

85= 3-77 „ -2-20, -1-60. 
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Prop. To find the difference in length, and also the dis- 
tance at the middle latitude, of two elliptic arcs of small 
ellipticity, lying in the plane of the meridian, and having 
their extremities in the same points; the latitudes of those points 
being known approximately, and the ellipses to which the arcs 
belong having their axes parallel. 

114. Let a and /3 be the co-ordinates to the centre of the 
ellipse, of which a, b, e are the semi-axes and ellipticity, 
measured from some fixed point near the centres of the two 
ellipses. The squares and products of a — b, e, a and /3 will 
be neglected. Let s be the length of the elliptic arc between 
the stations, l and l' the observed (or approximate) latitudes of 
the extremities, X and m the amplitude and middle latitude. 

First. We will find the length of the arc. Let c be the 
chord, r and 9, r and 9' the polar co-ordinates from the centre 
of the ellipse to the extremities of the arc. 

.*. c 8 = r*+ r' 8 — 2 rr' cos ( 9 — 0') = 2 rr' {1 —cos (0 — fl')} + (r— r')*, 

r = a (1 — esin*Z), r' = a (1 — e sin 8 ?'). 

Also tan 0 = (1 — 2e) tan l, 0 = l — e sin 2 ? ; 

.*. 9 — & = X — 2e sin X cos 2m ; 

.*. 1 — cos (0— 0') = 1 — cos X — 2e sin 8 X cos 2m 

= 2 sin 8 ^ X {1 — 2e (1 + cos X) cos 2m ) ; 

u 


.*. c 8 = 4 a 8 sin 8 ^ X {1 — 2e (1 + cos X) cos 2m- e (sin 8 1 + sin 8 ?')} 

A % 

= 4a 8 sin 8 ^ X [1 — e (1 + (2 + cos X) cos 2m}] ; 

A 

.'. sin \ X = ^ T 1 + \ e [l + (2 + cos X) cos 2m}"| ; 


P. A. 


*8 
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V = s “ ' Ya + 1 * I 1 + ( 2 + 008 X > C °" 2 ”> TShf 


= sin" 1 J e {1 + (2 + cos A) cos 2m} tan -X. 

Id It ** 

/ 1 \ 3 

Now s = a f 1 — g 6 J — 2 ae s ‘ n ^ 003 2w b *>y Art. 96 » 

= a (2— e) sin -1 + ae {1 + (2+cos X) cos 2m} tan | X 


— - ae sin X cos 2m 
2 


= (a+b) sin -1 ^ + (a - 5) j 1 + \ (1 - cos X) cos 2m 

Taking the variations, c being constant, 

l_£_ 

2 a 


tan - X. 

jS 


... , _. c a + b c8a 
Ss=($a + $&)sm — - 


V 4a s — c s 

+ (8a - 85) j 1 + ^ (1 - cos X) cos 2mj tan ^ X. 

The terms being small we may approximate ; 

/. Ss= (8a + 8b) ix — 2tan^X.Sa 

+ {8a- 8b) 1 1 + i (1 - cos X) cos 2mj tan |x 

= (8a+8i) Q X - tan ^ x) + ( 8a - 8b) ^ tan ^X(l-cosX) cos2m. 

115. Secondly. The distance between the arcs. 

The equation to the local ellipse is 

(»-«)* ,(y 
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•\ x*+y* o t r* = a’ + 2ax + 2/% — 2e (a* — a*) 

= a 3 + 2aa cos 0 + 2a/9 sin 0 — 2 a’e sin* 5 ; 

.*. r = a + a cos 0 + /3 sin 0 — ae sin* 0. 

+ ' 

Let R, C, C' be the values of r at the mid-latitude and at 
the extremities of the arc ; 


R = a + a cos m + $ sin m — (a — b) sin* m, 

G = a + a cos l + /9 sin l — (a — b) sin’ l, 

C — a + a cos V + /3 sin V — [a — b) sin’ l\ 

Multiply by 1, M, and N; add, and make the coefficients of 
a and /3 vanish ; 

.*. cos m +M cos l + N cos l' = 0, sin m+M sin l +N sin l' — 0; 

M = - 8 l n ^ ~ ^ | = sec \ A = N; 

sm {l — l) 2 2 

R+MC+NC' 

= a (1 + M+ N) — (a — b) (sin’ m + M sin’ l + N sin’ l') 

= a (1 4- 2 M) — i (a — b) {1 — cos 2 m + 2 M (1 — cos A cos 2m)} 

= i (a + i) (1 + 2 M) + ^ (a — b) (1 + 2Jf cos X) cos 2m 

A A 

= ^ (a + ft) ^1 — sec | X^ + ^ (a — i) ^1 — sec i X cos A^ cos 2m. 

i 

Taking the variations, the distance required, or 8R, 

= | (8a + 8b) f 1 - sec ^ X^ 

+ 1 (8a — 8b) ^1 - sec ^ A cos cos 2m. 
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116. In order to apply these expressions for Ss and SR to 
the arc in question, we must put 

X = 11*. 27'. 11", 2 m = 47°. 34'. 25" ; 

.*. Ss = 0*0003397 Sa - 0-0010097 Sb, 

SR = 0-0025536 Sa - 0-0075756 Sb. 

Substituting the values of Sa and Sb, Art. 113, 

Arc I. Arc II. Arc III. 

S3 = - 0-0000936, 0-0009543, 0-0001977 miles, 

SR = - 0-0006492, 0-0071414, 0-0089798 ... 

These quantities are so small as to he practically insensible : 
the largest value of Ss being l foot in an arc of 800 miles, and 
the largest value of SR being less than 12 yards. 

The result is, that the differences which are found to exist, 
between the observed amplitudes of arcs and the same ampli- 
tudes calculated geodetically, can in no respect be accounted 
for by supposing the arcs to be curved differently from the 
mean ellipse; because, as the above calculation shows, the 
ellipses may differ considerably in form without producing a 
sensible effect upon the length of the arc. This conclusion 
differs from that come to in the first edition of this work. In 
that edition (Art. 104) the distance of one extremity of the 
arc from the centre of the ellipse was taken to be the same in 
the local and the mean ellipse ; that is, those ellipses were 
supposed to be concentric, which they need not be, and are 
seen from the above investigation not to be. 

Prop. The differences between the astronomical and geo- 
detical amplitudes of an arc of meridian arise solely from local 
attraction, and are an accurate measure of the differences of 
local attraction at the extremities of the arc. 

/ 

117. The truth of this Proposition appears from the last 
Article. But we will establish it further by ascertaining how 
great a departure from the mean ellipse may exist without its 
producing even 1" of difference only in the amplitudes, as 
measured by the heavens and by the Earth. 
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By Art. 116, as 1°= 69*5 miles, 

0-00033978a - 0*001009786 = 1" = 0*0193056 ; 

8a — 56*8 + 2*97 8b miles. 

The values of 8a and 8b which satisfy this equation and 
make the sum of their squares least, are 8a = 5, 8b = — 17 
miles. But this variation in b is greater even than the whole 
compression of the pole, which is only 13 miles, to say nothing 
of the value of 8a in addition. It may be safely concluded, 
therefore, that no hypothesis regarding the curvature of the 
Indian Arc will account for the defect of 3"*791 and the excess 
of 5''*236 in the geodetic measure, which Colonel Everest 
found in the two arcs between Kalianpur and Damargida 
and between Kaliana and Kalianpur, when compared with 
the astronomical latitudes. 

There is no other possible cause, but local attraction affect- 
ing the plumb-line and level. These errors, therefore, become 
the accurate measure of the differences of the resultant local 
attraction, arising from causes visible and hidden, at the ex- 
tremities of the arc. The effect of the two visible causes, the 
mountain-mass and the ocean, taken together is very well 
represented (as already explained) by 3''*82 and 13"*11. To 
change these to — 3"*79 and 5"*24 (the values obtained by 
Colonel Everest from the comparison of the arc with the hea- 
vens) we must suppose some invisible cause, counteracting 
the effects of both the mountains and the ocean, and dimi- 
nishing their combined effect in these two arcs respectively 
by 3"*82 + 3"*79 or 7"*61 and 13"*ll-5"*24 or 7"*87. These 
quantities are nearly equal, and point to some cause existing 
in the crust beneath near the middle of the arc, that is, in 
the neighbourhood of Kalianpur. That even a slight excess 
of density through a large space around Kaliana is capable of 
producing such an effect we have shown in Art. 65. An 
endless variety of other hypotheses may be conceived to pro- 
duce this result, e. g. a deficiency of density beneath the 
mountains, accompanied by a corresponding deficiency south 
of Damargida towards Cape Comorin and the ocean. This 
double hypothesis is not, however, so simple as the single 
hypothesis above given. Whatever may be the facts of the 
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case, this is certain, that the difference of the local attractions 
in the meridian at Damargida and Kalianpur is 3"’79 south 
and the difference at Kalianpur and Kaliana is 5"*24 north. 

Prop. To explain what effect local attraction will have upon 
the mapping of a country. 

118. From what goes before it is clear that although the 
elliptic elements of the actual arc between two places may 
differ considerably from those of the mean ellipse, no sensible 
error will thence arise if we calculate the latitudes geodetically 
and with the elements of the mean ellipse. 

If however the latitudes are laid down in a map from ob- 
servations of the sun or stars they will be erroneous by the 
whole amount of deflection of the plumb-line by which the 
vertical and horizontal are determined. Thus in the case 
before us the deflections are as follows : — 


At Damargida, 

Kalianpur, 

Kalian a. 

By mountains 

6" -79 

12" -05 

27”*98 

„ ocean 

10 -44 

9 -00 

QO 

r— ! 

<x> 

„ hidden cause 

7 -61 

o 

© 

© 

- 7 *87 

Total deflections 

24 -84 

21 *05 

26 -29 


These angles converted into miles, at the rate of 1 # to 69‘5, 
or 51" '8 to 1 mile, are 0'48, 0‘41, 0'51 miles; by which quan- 
tities would the stations be wrongly placed on the map. The 
relative error is largest in the upper division of the arc, and 
in that case is not more than 1-1 Oth of a mile; but the positive 
error in each case is about half a mile. 

If, then, the principal places are all marked down geodeti- 
cally they will be correctly placed on the map, but if other 
places are filled in from observations of the sun or any other 
heavenly body they will be out of place by the whole of the 
error, viz. about half a mile. 

Prop. Geodesy furnishes no evidence, in proof or disproof, 
of the upheaval or depression of the Earth's surface as sug- 
gested by geological phenomena. 
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119. It will be observed, that the three arcs, which have 
been examined in Arts. 114 — 116, were compared, not with 
the mean ellipse itself, but with an ellipse equal in dimensions 
to the mean ellipse and with axes parallel (because the lati- 
tudes are measured in all the ellipses from the same or 
parallel lines). For this ellipse was so drawn as to pass 
through the extremities of the arc ; and we have no means of 
knowing that the mean ellipse passes through those two 
points. It may lie above them or below them. We have no 
means of ascertaining the position of the centre of the mean 
ellipse. The only way of doing this is to make a geodetic 
measurement of the whole of one meridian from pole to pole. 
Till this is done we have no evidence of any particular arc 
lying above or below the mean, i. e. of its having been 
elevated or depressed. The greatest geological changes of 
level, therefore, are perfectly consistent with all we know by 
geodesy of the surface of the Earth. 

120. In consequence of the inequalities of the Earth’s sur- 
face, levelling operations are all referred to the sea-level; 
that is, to that surface which the sea would form if allowed to 
percolate by canals through the continents. The sea is thus 
taken as the basis of our measurements ; and is assumed to 
have a spheroidal form. But it is possible that local disturb- 
ing forces, arising from attraction, may have the effect of 
crowding up the waters in the direction in which the forces 
act, so as sensibly to alter the sea-level from the spheroidal 
form. This we shall proceed to examine. 

Prop. To find the effect of a small horizontal disturbing 
force in changing the Level of the Sea. 


121. Let U be the disturbing force and du an element 
of the line u along which it acts. Then TJdu must be added 
to dV in the equation of fluid equilibrium of Art. 73. 



= V + 


ZD* f 

— r* (1 Udu = const, at the surface. 


Putting w' = m .E+e? and substituting for V (Art. 91); 
constant = f + (« “ f ) ^ (f - /**) + J f i 1 ~ ?*) + / Udu ‘ 
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When the small force U is neglected, a^-r = 1 + e . /**. 
Hence, neglecting small quantities of the second order, di- 
viding by E, multiplying by a, and transposing, 


* = constant + e./i* — Udu = 1 + e. /t * — ^ J 


Udu. 


1 dv 

Now - ^ is the tangent of the angle between r and the 

normal, = tan yjr suppose : and the angle through which the 
normal is thrown back by the force U 


Hence the element ds of the undisturbed meridian line on 
the surface of the sea is elevated, on the side towards which 
U acts, by the space 

=h u % da =-s Wu= j du - 

.*. whole elevation of the sea-level = - 

3 

integrated between the limits. 



ds . d^jr 


122. Ex. L. The Himmalayas attract places along the 
coast of Hindostan with a force varying nearly inversely as 
the distance from a line running E.S.E. and W.N.W. through 
a point in latitude 33° and longitude 77° 42', and equal to 
g tan 7" at 1020 miles distance: (see Phil. Trans. 1855, p. 91* 
94 ; also 1859, p. 793). Find the effect upon the sea-level 
between Cape Comorin and Karachi, which are about 1600 
and 775 miles from this line, arising from this cause. 

In this case U— — g tan 2" (1060 -r- u) u is the distance 
from the line. We may take the arc for the chord. There- 
fore rise of sea-level from this cause 

= 1020 tan Z" log, miles = 0-0346 x ~~~~ 

— 0-025 mile = 132 feet. 
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Ex. 2. As the distance from the line increases the force 
will vary more as the inverse square. Suppose that to the 
distance 1020 miles it varies as the inverse distance, and 
beyond that as the inverse square. For the first we must 
integrate as above : thus 

0-0346 log,, =0-0346 ^ = 0-0095 mile = 50 feet. 

775 0*434 

For the more southern part U=—g tan 7" (1020 -s- u) 2 , and 
the rise of the level 

=1020 tan 7 " f Wo346 x ^ = O'Ol 254 mile =66 feet. 

\1020 1600/ oO 

The sum of these is 116 feet, and is somewhat less than 
the result before obtained. We shall not be above the mark, 
therefore, in using the latter. 

Ex. 3. If u be the distance, in linear degrees, of the 
parallel of any place on the west coast of Hindostan from that 
of Cape Comorin, then the force acting towards the north at 
any point of that coast, arising from the deficiency of matter 
in the Ocean, may be approximately represented by the follow- 
ing formula (see Phil. Trans. 1859) : 

(0-000059556839 -0-000002836162m + 0-00000000407 2m") g. 

Hence at this place the sea-level is higher than at Cape 
Comorin, in consequence of this cause, by 

0-000059556839M - 0-000001418081m* + 0-000000001 357m 3 . 

Karachi is about 17° north of Cape Comorin. Hence, from 
this cause, the sea is higher at Karachi than at Cape Comorin 
0-00122 of a linear degree = 0-8489 mile = 448 feet. 

Ex. 4. Suppose an attracting force resides in Kalianpur, 
sufficient to produce a deflection 7"* 75 at 400 miles’ distance, 
and that the force varies inversely as the square of the dis- 
tance ; find its effect on the level between Cape Comorin and 
Karachi. 

p. a. 9 
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In this case the force is g tan 7" - 75 (400 -s- u)*, u being ex- 
pressed in miles ; 

rise in level = 160000 tan 7" • 76 = 21 feet, 

V«I V 

u 0 and u L being 1125 and 643 miles, the distances of Kalian- 
pur from Cape Comorin and Karachi. 

Taking tliQ sum of these three causes together, the increase 
in height of the sea- level at Karachi above that at Cape 
Comorin is 116 + 448 + 21 = 585 feet. There may be also 
other causes which may increase or decrease this result. But 
it serves to illustrate to what extent local attraction may have 
an effect upon the standard level to which all heights are 
referred. 


FINIS. 

I 


ERRATA. 

Page a, 1 . 2, for * read I 
„ 27, 1 . 4, insert the powers of c as in Art. 29 
„ 51, 1. 2, ab imo, for 12 read 21, for ^ read - sin 
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